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Abstract. We show that the regularity of the gravitational metric tensor 
in spherically symmetric spacetimes cannot be lifted from C"^^ to C^'^ by 
any C^'^ coordinate transformation in a neighborhood of a point of shock 
wave interaction in General Relativity, without forcing the determinant of the 
metric tensor to vanish at the point of interaction. This is in contrast to Israel's 
Theorem which states that such coordinate transformations always exist in a 
neighborhood of a point on a smooth single shock surface ^lOj . The results thus 
hiiply that points of shock wave interaction represent a new kind of spacetime 
singularity for perfect fluids, singularities that lie in physical spacetime, that 
can form from the evolution of smooth initial data, but at which the spacetime 
' ^ ' is not locally Minkowskian under any coordinate transformation. In particular, 

at such regularity singularities^ delta function sources in the second derivatives 
^ of the gravitational metric tensor exist in all coordinate systems, but due to 

j^JT) cancelation, the curvature tensor remains bounded (in L°°). 

' We announced the main results of this paper in [151 , where we also sketched 

the main steps in some of the proofs. In this paper, those results together with 
^ their complete proofs and some additional propositions are presented. 
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1. Introduction 

In 1915, seeking to conciliate Newtonian gravity with Special Relativity, Al- 
bert Einstein formulated the theory of General Relativity [3]. At its heart lies the 
(gravitational) metric tensor describing the geometry of spacetime and thus the 
gravitational field. A novel feature of General Relativity (GR) was the predic- 
tion of spacetime singularities, those are points where the gravitational metric 
suffers a lack of regularity. For instance, the singularity at the center of the 
Schwarzschild metric or at its Schwarzschild radius, where the metric tensor fails 
to be bounded. The first one is an example of a non-removable singularity which 
persist in every coordinate system, those singularities are usually characterized 
by a blow-up in the scalar curvature and they lie outside of physical spacetime. 
The apparent singularity at the Schwarzschild radius is an example of a remov- 
able singularity, i.e., there exist coordinates in which the metric is regular enough 
to be non-singular. A metric is non-singular if its components, their first and sec- 
ond derivatives and its curvature are bounded and if it is locally inertialj^that is, 
around any point p exist coordinates in which the gravitational metric at p is the 
Minkowski metric up to second order corrections. (The physical meaning of the 
metric being locally inertial is that an observer in freefall through a gravitational 
field should observe the physics of special relativity, except for the second order 
acceleration effects due to the gravitational field.) However, the metric tensor is 
governed by the Einstein equations, which are a system of PDE's, so that the 
Einstein equations by themselves determine the smoothness of the gravitational 
metric tensor by the evolution they impose. Thus the condition on spacetime 
that it be non-singular cannot be assumed at the start, but must be determined 
by regularity theorems for the Einstein equations. 

Whenever the sources of matter and energy are modeled by a perfect fluid, 
this issue becomes all the more interesting and intriguing. Then the Einstein 
equations imply the GR compressible Euler equations through the Bianchi iden- 
tities, and the compressible Euler equations create shock waves out of smooth 
initial data whenever the flow is sufficiently compressive [21 [161 [Ej- At a shock 
wave, the fluid density, pressure, velocity, and hence T are discontinuous, so 
the Einstein equations imply the curvature G must also become discontinuous 
at shocks. But discontinuous curvature by itself is not inconsistent with the 
assumption that spacetime be non-singular. For example, if the gravitational 
metric tensor were C^'^, (differentiable with first derivatives being Lipschitz con- 
tinuous, C°'^), then second derivatives of the metric are at worst discontinuous, 
the curvature is bounded (in L°°) and the metric has enough smoothness to be 
locally inertial, |18]. Furthermore, Israel's theorem asserts that a metric C°'^ 
regular across a smooth single shock surface, is lifted to G^'^ by the G^'^ coor- 
dinate map to Gaussian normal coordinates, and this is again smooth enough 
for spacetime to be non-singular at each point. In [7j, Groah and Temple set 
out a framework in which to address these issues rigorously by providing the 
first general existence theory for spherically symmetric shock wave solutions of 
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the Einstein- Euler equations allowing for arbitrary numbers of interacting shock 
waves of arbitrary strength. In coordinates where their analysis is feasible, Stan- 
dard Schwarzschild Coordinates (SSC), (a general spherically symmetric metric 
can generically be transformed to SSC, c.f. [25]), the gravitational metric is only 
C°'^ at shock waves, and it has remained an open problem as to whether the 
weak solutions constructed by Groah and Temple could be smoothed to C^'^ by 
coordinate transformation, like the single shock surfaces addressed by Israel. 

The resolution of the open problem of Groah and Temple is achieved in [15] 
and the present paper by proving that there do not exist C^'^ coordinate trans- 
formations that can lift the regularity of a gravitational metric tensor from C°'^ 
to C^'^ at a point of shock wave interaction in a spherically symmetric spacetime, 
without forcing the determinant of the metric tensor to vanish at the point of 
interaction. Consequently, in contrast to Israel's theorem for single shock sur- 
faces, shock wave solutions cannot be continued as C^'^ strong solutions of the 
Einstein equations beyond the first point of shock wave interaction. The results 
imply that points of shock wave interaction represent a new kind of singularity 
in General Relativity that can form from the evolution of smooth initial data, 
that lies within physical spacetime, but at which second order metric derivatives 
are distributional and spacetime is not locally inertial under any C^'^ coordinate 
transformation. Due to cancelation, the Einstein tensor is sup-norm hounded 
and free of delta function sources at points of shock wave interaction. This re- 
sult contrasts the common assumption about the metric being C^'^ regular, for 
example, this is assumed in the singularity theorems of Hawking and Penrose, 
[9]. In this paper we present the complete proofs of the results announced in 
|15j . where many proofs were omitted. 

To state the main result precisely, let g^j^y denote a spherically symmetric space- 
time metric in SSC, that is, the metric takes the form 

rfs^ = g^.udx'^dx" = -A{t, r)dt'^ + B{t, r)dr^ + r^dQ^, (1.1) 

where either t or r can be taken to be timelike, and dfl"^ = d'd'^ + siv?['d)dip'^ is 
the line element on the unit 2-sphere, c.f. [19]. In Section [s] we make precise the 
definition of a point of regular shock wave interaction in SSC. Essentially, this is 
a point where two shock waves enter or leave the point p at distinct speeds, such 
that the metric is Lipschitz continuous across each shock, the Rankine Hugoniot 
(RH) jump conditions hold across the shocks, and the SSC Einstein equations 
hold strongly away from the shocks. The main result of this paper is the following 
theorem, c.f. [15] : 

Theorem 1.1. Assume p is a point of regular shock wave interaction in SSC, in 



the sense of Definition 3.1, for the SSC metric gfj_u- Then there does not exist a 
C^'^ regular coordinate transformation, defined in a neighborhood of p, such that 
the metric components are functions of the new coordinates and such that the 
metric has a nonzero determinant at p. 



The proof of Theorem 1.1 is constructive in the sense that we characterize the 



Jacobians of (t, r) coordinate transformations that could smooth the components 
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of the gravitational metric in a deleted neighborhood of a point p of regular shock 
wave interaction, and then prove that any such Jacobian must have a vanish- 
ing determinant at p itself. We then extend the result to C^'^ transformations 
allowing for changes of angular variables as well. 

We emphasize, when discontinuities in the fluid are present, C^'^ coordinate 
transformations constitute the atlas of transformations capable of lifting the 
regularity of the metric one order (across a single shock surface), while still 
preserving the weak formulation of the Einstein equations, c.f. |T8]. For this 
and other reasons, c.f. Section [lT| is the atlas of C^'^ coordinate transformation 
generic for addressing shock wave interactions in GR. 

Our assumptions in Theorem 1 1 . 1 1 apply to the upper half {t > 0) and the lower 
half (if: < 0) of a shock wave interaction (at t = 0) separately, general enough to 
include the case of two timelike (or spacelike) interacting shock waves of opposite 
families that cross at the point p, but also general enough to include the cases 
of two outgoing shock waves created by the focusing of compressive rarefaction 
waves, or two incoming shock waves of the same family that interact at p to create 
an outgoing shock wave of the same family and an outgoing rarefaction wave of 
the opposite family, c.f. [T7]. In particular, our framework is general enough to 
incorporate the shock wave interaction which was numerically simulated in [24j . 

Even though our research was motivated by shock wave solutions of the Einstein- 
Euler equations, we do not assume a perfect fluid as the matter source in The- 
orem |l.l[ so that our main result is applicable to other matter models as well, 
provided T^^ is bounded and satisfies the RH conditions, \T^^^]n^ = 0, on each 
of the shock curves. 

Historically, the issue of the smoothness of the gravitational metric tensor 
across interfaces began with the matching of the interior Schwarzschild solution 
to the vacuum across an interface, followed by the celebrated work of Oppen- 
heimer and Snyder who gave the first dynamical model of gravitational collapse 
by matching a pressureless fluid sphere to the Schwarzschild vacuum spacetime 
across a dynamical interface [13]. In [18], Smoller and Temple extended the 
Oppenheimer-Snyder model to nonzero pressure by matching the Friedmann 
metric to a static fluid sphere across a shock wave interface that modeled a 
blast wave in GR. In his celebrated 1966 paper |10J, Israel gave the definitive 
conditions for regular matching of gravitational metrics at smooth interfaces, by 
showing that if the second fundamental form is continuous across a single smooth 
interface, then the RH conditions also hold, and Gaussian normal coordinates 
provide a locally inertial coordinate system at each point on the surface. In [7| 
Groah and Temple addressed these issues rigorously in the first general existence 
theory for shock wave solutions of the Einstein-Euler equations in spherically 
symmetric spacetimes. The results were extended to Gowdy spacetime in fT]. 

Although points of shock wave interaction are straightforward to construct 
for the relativistic compressible Euler equations in flat spacetime, we know of 
no rigorous construction of a point of regular shock wave interaction in GR. 
However, all evidence indicates points of shock wave interaction to exist, have 
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the structure we assume in SSC, and cannot be avoided in solutions consisting of, 
say, an outgoing spherical shock wave (the blast wave of an explosion) evolving 
inside an incoming spherical shock wave (the leading edge of an implosion). 
Namely, the existence theory of Temple and Groah [7j lends strong support to this 
claim, establishing existence of weak solutions of the Einstein- Euler equations in 
spherically symmetric spacetimes. The theory applies to arbitrary numbers of 
initial shock waves of arbitrary strength, existence is established beyond the point 
of shock wave interaction, and the regularity assumptions of our theorem are 
within the regularity class to which the Groah- Temple theory applies. Moreover, 
the recent work of Vogler and Temple gives a numerical simulation in which two 
shock waves emerge from a point of interaction where two compression waves 
focus into a discontinuity in density and velocity, and the numerics demonstrate 
that the structure of the emerging shock waves meet the assumptions of our 
theorem. 

It is instructive at this point to clarify the difference between the essential 
(jo,i singularities in the metric at points of shock wave interaction, and the 
essential singularities at surface layers like the "thin shells" introduced in 
Israel's illuminating paper [IDj. (See also |B].) On surface layers, the delta 
function sources in T are the cause of the essential C^'^ singularity in the metric g, 
because second derivatives of g must have distributional sources and consequently 
g cannot be C^'^ in any regular coordinate system. For shock wave solutions of 
G = kT the issue is more delicate because T is sup-norm bounded, so that the 
constraint of G having delta function sources is removed and, at first sight, there 
is no clear obstacle to the existence of coordinate systems that smooth the metric 
to C^'^. Israel's theorem confirms there is no obstacle to C^'^ smoothness in the 
special case of single shock surfaces, but the methods in [7] are only sufficient 
to prove existence of solutions in C*^'^, and the question as to whether there 
is an obstacle for more complicated solutions with interactions has remained 
unresolved until now. 

The argument in [15] and the present paper resolves this issue by proving 
that at points of shock wave interaction, the Einstein-Euler equations in SSC 
generate an essential C°'^ singularity in the metric that cannot be smoothed 
to C^'^ by coordinate transformation, even though there are no delta function 
sources present in T or G to explain the C'''^ singularity in the metric g. 

We conclude that points of shock wave interaction are a new kind of regularity 
singularity in the gravitational field that are not generated by delta function 
sources in T, that can form from the evolution of smooth initial data, that 
correctly refiect the physics of the equations even though the spacetime is not 
locally Minkowskian under any local C^'^ coordinate transformation, and where 
the metric tensor does not have sufficient regularity to satisfy the Einstein-Euler 
equations strongly in any coordinate system of the C^'^ atlas. At such singular- 
ities, delta function sources in the second derivatives of the gravitational metric 
tensor exist in all coordinate systems of the C^'^ atlas, but due to cancelation, 
the Riemann curvature tensor remains sup-norm bounded. 
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In Section [3} we set up our basic framework and define the notion of a point 
of regular shock wave interaction in SSC We make precise what we mean by a 
function (or a metric) to be C°'^ across a hypersurface and derive a canonical 
form for such functions in Section |4j In Section [5} we show that the property 
of the metric tensor being C°'^ across a shock surface is covariant under 
coordinate transformations, but not under C^'^ transformations. We use this 
lack of covariance to derive conditions on the Jacobians of general C7^'^ coordinate 
transformations (in principal) necessary and sufficient to lift the regularity of a 
metric tensor from C°'^ to at points on a shock surface. This condition enables 
us to represent all such Jacobians in terms of the canonical form introduced in 
Section |4| unique up to addition of an function. 

In Section [6} we give a new constructive proof of Israel's theorem for spher- 
ically symmetric spacetimes, by showing that the Jacobians expressed in our 
canonical form do smooth the gravitational metric to C^'^ at points on a single 
shock surface. The essential difficulty is to prove that the "gauge" freedom to 
add an arbitrary C^-function to our canonical form is sufficient to satisfy the in- 
tegrability condition, which is required to integrate the Jacobian to coordinates. 
This is possible within the gauge freedom if and only if the so-called Rankine 
Hugoniot conditions and the Einstein equations hold at the shock interface. 

The main step towards Theorem 1.1 is achieved in Section [7] where we prove 
that at a point of regular shock wave interaction in SSC there exist no coordinate 
transformations in the [t, r)-plane, (that is, transformations that keep the angu- 
lar part fixed), that lift the metric regularity to C^. The central method herein 
is that the gauge freedom in our canonical forms is insufficient to satisfy the 
integrability condition on the Jacobians, without forcing the determinant of the 
Jacobian to vanish at the point of interaction. In Section [8] we extend this result 
to the full atlas of spacetime, thereby proving Theorem |1.1[ We prove the loss of 
locally inertial frames in Section |9] and show that the Riemann curvature tensor 



is bounded at points of shock wave interaction in Section 10 



2. Preliminaries 

In this section we discuss the framework of General Relativity and introduce 
the coupled Einstein-Euler equations together with shock waves as solutions of 
the latter. In addition, we give a short review of Israel's result [10', T!8], introduce 
spherically symmetric spacetimes and discuss the Einstein equations in Standard 
Schwarzschild Coordinates regarding the metric regularity. 

Spacetime is a four dimensional Lorentz manifold, that is, a manifold endowed 
with a metric tensor of signature (—1, 1, 1, 1), the so-called Lorentz metric. Re- 
quiring spacetime to be a manifold reflects Einstein's original insight of general 
covariance [3j, that is, all physical equations must be formulated as tensor equa- 



tions (c.f. (2.2)). The signature of the metric induces a notion of causality 
P[261E5]. 

For this paper, an essential feature of a manifold is that around each point 
p E M there exists an open neighborhood JVp (called a patch) together with 
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a homeomorphism ^) : A/'p — !■ M^, which defines coordinates on 

its image in M^. Mp together with the homeomorphism x are called a chart. 
The collection of all such charts (covering the manifold) is called an atlas. If 
the intersection of two coordinate patches is nonempty, x o defines a map- 
ping from an open set in to another one (often referred to as a coordinate 
transformation). Given that all coordinate transformations in the atlas are 
differentiable, the manifold M is called a C'^-manifold and its atlas a C'^-atlas. In 
General Relativity the manifold is usually assumed to be C^, however, in order 
to address shock wave solutions of the Einstein Euler equations it is crucial to 
consider C^'^-manifolds. In fact, lowering the regularity to C^'^ is the crucial step 
allowing for a smoothing of the metric in the presence of a single shock wave, 
(c.f. [sjand Theorem [gJ. 

Throughout this paper we use the Einstein summation convention, that is, we 
always sum over repeated upper and lower indices, e.g. v^w^ = v^wq + ... + v^ws 
for /i, G {0,...,3}. Moreover, we use the type of index to indicate in which 
coordinates a tensor is expressed in, for instance, T^^, denotes a (l,l)-tensor 
in coordinates x'^ and T"^ denotes the same tensor in (different) coordinates 
x". Under a change of coordinates, tensors transform via contraction with the 
Jacohian of the coordinate transformation, which is given by 

= — (2 1) 

dx- ■ ^ ^ ^ 

For example, a vector transforms as = J^f", a one-form as = J^v^ and a 



(1, l)-tensor as T'^j, = J^J^T^p, where denotes the inverse of (2.1), c.f. 
for some basics on tensors. The metric tensor transforms according to 

9^.u = r.J^ap , (2.2) 

which is crucial for the methods in this paper. We denote with g^^ the inverse of 
the metric, defined via 

g^'^ga. = 5^, (2.3) 

in terms of the Kronecker-symbol 6^. By convention, we raise and lower tensor 
indices with the metric, for example T^,^ = g^^rT'^'^. 



In (2.1), the Jacobian is defined for a given change of coordinates. Vice versa, 
for a set of given functions to be a Jacobian of a coordinate transformation 
it is necessary and sufficient to satisfy (see appendix |A] for details), 

= (2.4) 
det(J„^) ^ 0, (2.5) 

where = denotes partial differentiation with respect to coordinates x". 



(2.4) ensures to be integrable to coordinates x^, while (2.5) ensures the 



invertibility of the coordinate functions and that non-singular tensors remain 



non-singular under change of coordinates. We henceforth refer to (2.4) as the 
integrability condition. 

We now discuss the Einstein equations that govern the gravitational field 
through coupling the spacetime curvature to the energy- and matter-content. 
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The Einstein equations read jSj 

Qi^u = ^^T^'^ (2.6) 

(in units where c = 1), where k := Svr^ incorporates Newton's gravitational 
constant Q into the equation and 

Qt^v ^ j^t^v _ ^J^g^,u ^ J^g^,u (2.7) 

is the Einstein tensor, A denotes the cosmological constant. (Our method ap- 
phes regardless the choice of A G M, since the term Ag^^^ gives a continuous 
contribution to the Einstein equations ( 2.6[ ).) The Ricci tensor R^i, is defined 
to be the trace of the Riemann tensor = W^^,^ and the trace of the Ricci 
tensor R = is called the scalar curvature. The metric tensor g^j^^ enters the 
Einstein equations through the Christoffel symbols 

1 

2' 

since the Riemann curvature tensor is given by 

DM = _ F/^ _|_ _ F/^ 

^vpu I'p,^ i"^,P Act up \p ua ' 

Taking into account the symmetry of the metric tensor g^^iy = g,yn, (and the Ricci 
tensor R^i, = R^p), the Einstein equations form a set of 10 second order differen- 
tial equations on g^j,^. For the Einstein tensor to be defined in a strong (almost 
everywhere) sense a C^'^ metric regularity is necessary, however, at the level of 
shock waves only Lipschitz continuity is guaranteed and one must introduce the 
Einstein tensor in a weak (distributional) sense, c.f. |18j . 

By construction, due to the Bianchi identities of the Riemann curvature tensor, 
the Einstein tensor is divergence free 

G"'',^ = 0, (2.9) 



^up o^'^ i9ucr,p + 9pcr,i> 9vp,a) i (2-^ 



where the semicolon in (2.9) denotes covariant differentiation, that is 

<. = <. + r^X (2.10) 



for a vector field v^. Through (2.9), the Einstein equations ensure conservation 
of energy in the matter source, that is, 

T^^i. = 0, (2.11) 

which was one of the guiding principles Einstein followed in the construction 
of the Einstein tensor [3]. In the case of a perfect fiuid, that is, the energy 
momentum tensor is given by 

T^"' = {p + p)u^'u'' + pg^"", (2.12) 



(2.11) are the general relativistic Euler equations, with p being the density, p 
the pressure and the tangent vector of the fiuid fiow normalized such that 
u^Uy = — 1 ,(see for example [20] )• In a locally inertial frame around a point p, 
that is, in coordinates where the metric satisfies (c.f. Definition 9.1) 

9pv{p) = Vpu 
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and gf,uAp) = V/i, z/,0- G {0,...,3}, (2.13) 
with r]^^ the Minkowski metric, ( |2.11 ) reduce to the special relativistic Euler 



equations, 



at the point p. (2.6) together with (2.11), are the coupled Einstein Euler equa- 



tions. Having unknowns g^^, P and u''', the coupled Einstein Euler equations 
are under-determined, however, prescribing an (barotropic) equation of state, 
P = P{p)y the system closes. 



The Euler equations (2.11) are a system of conservation laws, where shock 
waves form from smooth initial whenever the flow is sufficiently compressive, [21 
[121 [17] . This makes the study of shock waves inevitable for perfect fluid sources, 
arguing the result of this paper being fundamental to General Relativity. 

Being discontinuous, shock wave solutions satisfy the conservation law only 
in a weak sense. To obtain the weak formalism of a PDE, first multiply the 
equation with a smooth ("test") function, then integrate the resulting equation 
and finally apply the divergence theorem to shift all derivatives onto the test 
function, c.f. [IH]. The weak form of the Euler equations reads 

T^V,.c?/^M = 0, (2.14) 



M 



where (p G C^{M) is a test function and dpM is the spacetime volume element. 
Across the surface of discontinuity S (the so-called shock surface) the solution 
satisfies the Rankine Hugoniot conditions (RH conditions), 

[T^"']N^ = 0, (2.15) 

where N'^ is normal to the hypersurface S and [u] := — ur for Ul and Ur 
denoting the left and right limit (to S) respectively, [u] is referred to as the 
"jump" in u across S and provides a measure for the strength of the shock. 
In fact, suppose T^*^ is a strong solution everywhere away from the surface of 
discontinuity, then T'^'^ is a weak solution in a whole neighborhood of the surface 



if and only if the RH conditions (2.15) hold everywhere on the surface. Thus, 



one can bypass the weak formalism and instead work with the RH conditions 
and the strong solution away from the shock surface, as we do in this paper. 
Again, we point out that we do not specifically use perfect fiuid sources in our 



method, it suffices to assume that T^'^ satisfies the RH condition (2.15) on S 



and is continuous away from some S, (c.f. Definition 3.1). 

For the Einstein tensor to be defined in a strong sense almost everywhere a C^'^ 
metric-regularity is required. However, at the level of shock waves only Lipschitz 
continuity of the matric is guaranteed, as we discuss below (see also p^J). For 
this low regularity objects like the Einstein or Riemann curvature tensor can 
only be introduced in a weak (distributional) sense. Nevertheless, in his 1966 
paper [10] Israel proved that, at a single smooth shock surface, there exists 
coordinates where metric is C^'^ regular if and only if the energy momentum 
tensor is bounded almost everywhere. If any of the above holds, the Rankine 
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Hugoniot jump conditions are satisfied everywhere on the shock surface. Another 
part of Israel's precise result is that the gravitational metric is smoothed to C^'^ 
in Gaussian Normal Coordinates if and only if the second fundamental form of 
the metric is continuous across the surface. The latter is an invariant condition 
meaningful for metrics Lipschitz continuous across a hypersurface, and is often 
referred to in the literature as the junction condition, c.f. [25] . 

At the heart of Israel's method lies the ad-hoc choice of Gaussian Normal 
Coordinates with respect to the shock surface E, that is, we first arrange by a 
smooth coordinate transformation that locally S = {p G M : = 0} and 

then exchange the n — th coordinate function by arc-length of geodesies normal 
to S. In more detail, Gaussian Normal Coordinates are the mapping assigning 
a point p G M (sufficiently close to S) to a point x"{p) in M" as follows: 

x"ip) = {s,x--\q),...,x\q)), (2.16) 

where s is the arc-length parameter of a geodesic curve 7 starting at the point 
g G S in the direction normal to S, with 7(5) = p, and x"(g) = x*(g) for 
all a = i G {l,...,n — 1}. Computing now the Einstein tensor in coordinates 



(2.16), one finds that each component of the resulting Einstein tensor contains 
only a single second order normal derivative, gaf5,nn, while all other terms in the 
Einstein equations are in L°° and thus ga/5,nn £ L°° as well. Since all other second 
order metric derivatives are bounded by assumption we conclud^^fa^ G C^'^ 
In Section [6] we give a new constructive proof of Israel's result, based on the 
method we introduce in Section [5l 

In the following we introduce and discuss spherically symmetric spacetimes. 
Many spacetimes of fundamental interest to General Relativity are spherically 
symmetric, for example, the Schwarzschild, the Oppenheimer Volkoff or the 
Friedmann Robertson Walker spacetimes, [H]. A spherically symmetric space- 
time is a Lorentz manifold, allowing for three independent spacelike Killing vec- 
tor field, such that the subspaces parameterized by the fiow of the Killing vectors 
have a positive constant curvature [25] . We refer to those subspaces as the spaces 
of symmetry, which in suitable coordinates are given by a family of two spheres 
of smoothly varying radii. A vector field is called a Killing vector if it satisfies 
Killing's equations 

X^., + X,.^ = 0. (2.17) 

Killing's equation ensures that the fiow of a solution X'^ is an isometry of space- 
time, since a vector field X^ solves ( 2.17[ ) if and only if the Lie derivative of the 



metric in direction of X'^ vanishes [9] [261 ES] , 

Lx9 = 0. (2.18) 

In a spherically symmetric spacetime, assuming that one of the spaces of sym- 
metry has constant scalar curvature K = 1, one can always introduce coordinates 



^Here we use that Lipschitz continuity of a function is equivalent for it to be in the Sobolev- 
space W^'°°, containing all functions with first (weak) derivatives in L°° . Both imply the 
function to be differentiable almost everywhere [4j. 
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and <y9 such that the metric takes on the form 

ds^ = -Adt^ + Bdr"^ + 2Edtdr + Cd^^ 



(2.19) 



where 



dn^ 



dip'^ + sin^ (^9) d^^ 



is the hne element on the two-sphere and the metric coefficients A, B, C and 
E only depend on t and r, [25]. This simplifies the metric significantly, as its 
original ten free components reduce to four and the ten Einstein equations reduce 



to four independent ones accordingly. The metric structure (2.19) is preserved 



under coordinate transformations in the (t, r)-plane, that is, transformations 
which keep the angular variables d and 99 fixed. Under generic conditions, one 
can simplify the metric further by introducing a new "radial" variable r' := 
\/~C and removing the off-diagonal element through an appropriate coordinate 
transformation in the (t,r')-plane denoting the resulting coordinates again 
by t and r the new metric reads 

(2.20) 



ds^ = -Adt^ + Bdr"^ + r^dfi^ 



Coordinates where the metric is given by (2.20 ) are called Standard Schwarzschild 



Coordinates (SSC). In SSC, the metric has only two free components and the 
Einstein equations simplify significantly: 

B - 1 



Br + B- 



with 



y4.^ 
Bu 

BAtBt 
2AB 



1 



Bt 
B 



A- 

r 

A„. + $ 



-kAB^T'^^ 

kAB\T^^ 

-2KABr^T^^ 



(2.21) 
(2.22) 
(2.23) 
(2.24) 



Bl 
2B 



— + 



rB 



A"^ 



2B 



The first three Einstein equations in SSC are central in the method we develop 



and in the proof of Theorem 7.1 



Finally, we address the issue of the metric regularity in the presence of shock 
waves. A shock wave is a (weak) solution T^'^ of the Euler equations (2.11), 



discontinuous across a timelike hypersurface S and smooth away from S, where 



the discontinuity satisfies the RH conditions (2.15) across S. Now, from the first 



three Einstein equations in SSC, (2.21)-(2.23), it is straightforward to read off 



that the metric cannot be any smoother than Lipschitz continuous if the matter 
source T'^'^ G L°° is discontinuous. In this paper we henceforth assume that the 
gravitational metric in SSC is Lipschitz continuous]^ providing us a consistent 
framework to address shock waves in General Relativity, agreeing with various 
examples of solutions to the coupled Einstein Euler equations, for example, the 



It is not clear to us if shock wave solutions with a lower metric regularity do exist. 
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solutions in [7] and [TH]. Moreover, Lipschitz continuity arises naturally in the 
problem of matching two spacetimes across an interface pUj . 

As mentioned above, Israel proved the remarkable result that whenever a 
metric is Lipschitz continuous across a smooth single shock surface S and has an 
almost everywhere bounded Einstein tensor, then there always exists a coordinate 
transformation defined in a neighborhood of E smoothing the components of the 
gravitational metric to C^'^. In spherically symmetric spacetimes, SmoUer and 
Temple showed in [18] that in a neighborhood of a single radial shock surface 
there exists coordinates where the metric is C^'^ regular if and only if the RH 
conditions (2.15) hold everywhere on the surface. Thus, for example, the metric 
across single radial shock surfaces can be no smoother than Lipschitz continuous 
in SSC, but can be smoothed to C^'^ by coordinate transformations. However, 
it has remained an open problem whether or not such a theorem applies to the 
more complicated C'^'^-solutions containing shock wave interactions, for example, 
the SSC solutions proven to exist by Groah and Temple [7j. Our purpose here 
is to show that such solutions cannot be smoothed to in a neighborhood of a 
point of regular shock wave interaction in SSC, a notion we now make precise. 



3. A Point of Regular Shock Wave Interaction in SSC 

In this section we set up our basic framework and give the definition of a point 
of regular shock wave interaction in SSC, that is, a point p where two radial shock 
waves enter or leave the point p with distinct speeds. We henceforth restrict 
attention to radial shock waves, that is, shock surfaces E that can (locally) be 
parameterized by 

E(t,^9,y.) = (t,x(t),^9,</.), (3.1) 

and across which the stress-energy-momentum tensor T is discontinuous. Said 
differently, for a fixed time t the shock surface is a two sphere (of symmetry) with 
radius x(t). Our subsequent methods apply to spacelike and timelike surfaces 
alike, (inside or outside a black hole, c.f. [T^), since we did not need to specify 
the signs of the metric coefficient, but without loss of generality and for ease of 
notation we henceforth restrict to timelike surfaces. 

For radial hypersurfaces in SSC, the angular variables play a passive role, and 
the essential issue regarding smoothing the metric components lies within the 
atlas of C^'^ coordinate transformations acting on the (t,r)-plane, i.e., angular 
coordinates are kept fixed. (In fact, for the proof of Theorem 6.1 and 7.1 we 
consider only the [t, r)-plane and standard methods suffice to extend the result 
to the full atlas, leading to Theorem |1.1[ ) Therefore it suffices to work with the 
so-called shock curve 7, that is, the shock surface E restricted to the {t, r)-plane, 

7(t) = (t,x(t)), (3.2) 

with normal 1-form 

= (x, -1). (3.3) 
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For radial shock surfaces (3.1) in SSC, the RH conditions (2.15) take the simph- 
fied form 

[T°°]x = [roi], (3.4) 

[T'^]x = [T'']. (3.5) 

In Section [7] we prove the main result of this paper by establishing that SSC 
metrics Lipschitz continuous across two intersecting shock curves, cannot be 
smoothed to C^'^ by a C^'^ coordinate transformation. To establish our basic 
framework, suppose two timelike shock surfaces Sj are parameterized in SSC by 

E,(t,^,0) = (t,a;,(t),e,0), (3.6) 

for i = 1,2, 3, 4, where Si and S2 are defined for t < and S3 and S4 are defined 
for if: > Oj^ described in the (t,r)-plane by, 

7,(t) = (t,x,(t)), (3.7) 

with normal 1-forms 

(n,), = (Xi,-1). (3.8) 
Denoting with the jump across the i — th shock curve the RH conditions read, 



in correspondence to (3.4)-(3.5) 



[T°°]^i;, = [T°i]^, (3.9) 
= [T^% (3.10) 



For the proof of our main result (Theorem 7.1) it suffices to restrict attention 
to the lower (t < 0) or upper (t > 0) part of a shock wave interaction that occurs 
at t = 0. That is, it suffices to impose conditions on either the lower or upper 
half plane (c.f. Figure 1) 

= {(t,r) : t < 0}, 

or 

Rl = {it,r):t>0}, 
respectively, whichever half plane contains two shock waves that intersect at p 
with distinct speeds. Thus, without loss of generality, let t < and let 'ji(t) = 
{t,Xi{t)), {i = 1,2), be two shock curves in the lower (t,r)-plane intersecting in 
a point (0,ro), for tq > 0, that is, 

a;i(0) = ro = X2(0). (3.11) 

We are now prepared to give the definition of what we call a point of regular 
shock wave interaction in SSC. By this we mean a poinlj^ p where two shock 
waves collide with distinct speeds, such that the metric is smooth away from the 
shock curves and Lipschitz continuous across each shock curve, allowing for a 
discontinuous T^^" and the RH condition to hold. Recall, we assume without loss 
of generality a lower shock wave interaction in M^. 



^Shock waves typically change their speeds discontinuously at the point of interaction. 
^The intersection of the shock surfaces is a two sphere, but abusing language we refer to it 
as a point, consistent with us suppressing the (trivial) angular dependence in our method. 
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Figure 1. Example of intersecting shock curves. 



Definition 3.1. Let Vq > 0, and let g^i, he an SSC metric in C°'^ (^Af H R 

where Af G M."^ is a neighborhood of a point p = (0, ro) of intersection of two 
timelike shock curves ^i{t) = {t,Xi{t)) G M^, t G (— e,0). Assume the shock 
speeds xAO) = limxj(t) exist and are distinct, Xi(0) ^ X2(0), and let J\ denote 

the neighborhood consisting of all points in M H not in the closure of the two 
intersecting curves '~fi{t). Then we say that p is a point o/ regular shock wave 
interaction in SSC if: 



(i) The pair {g,T) is a strong solution of the SSC Einstein equations (2.21) 



( [2:241 ) «^ with T^"' G C^iAf) and g^^ G C^Af). 
(a) The limits ofT^'^ and of metric derivatives g^v^a exist on both sides of each 

shock curve ■jiit) for all G (— e, 0). 
(Hi) The jumps in the metric derivatives [gfj.u,o-]i{t) are function with respect 

to t for 2 = 1,2 and for t G (— e, 0) . 
(iv) The limits 

lim[^^^,<,]i(t) = [gf,^,a]i{0) 

exist for i = 1,2. 

(v) The metric g is continuous across each shock curve 'yi{t) separately, but no 
better than Lipschitz continuous in the sense that, for each i there exists 
fj,, u such that 

[g^.vA^{n^Y ^ 



denotes the closure of Mj. . 
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at each point on t G (— e, 0) and 

\m^g^y^„]i{niY ^ 0. 

(vi) The stress tensor T is hounded on M r\ M?_ and satisfies the RH conditions 

[T'^Uui)^ = 

at each point on '~fi{t), i = 1,2, t E (— e,0), and the limits of these jumps 
exist up to p as t — )■ 0. 



The structure assumed in Definition 3.1 reflects the regularity of shock wave 



solutions of the coupled Einstein-Euler equations consistent with the theory in 
[7] and conflrmed by the numerical simulation in [23] . If one wants to include the 
structure of single general relativistic shock waves, for instance, as in fl8{ [T9| l20t 



we expect Deflnition 3.1 to be the most natural way of introducing 
shock wave interactions. In more detail, (i)-(iii) and (vi) are straightforward 
generalizations from either the flat case or the case of a single relativistic shock 
wave. Note, we do not require T^^^ to match up continuously at p. Condition (iv) 
prohibits blow-up phenomena in the matter sources and condition (v) ensures 
the shock to persist up to the interaction. 

Our method of proof works for timelike and spacelike shock surfaces alike, but 
breaks down for null hypersurfaces. Also, our method fails in the unphysical case 



x^{0) = X2{0), due to (7.19). 



Note, we do not explicitly assume a perfect fluid source in Deflnition 3.1 
only an energy momentum tensor T^^" that satisfles the Rankine Hugoniot jump 
condition. Furthermore, the structure we deflned is general enough to include 
all shock wave interaction described in [1^. In more detail, it includes the case 
of two interacting shock waves of opposite families that cross at the point p, but 
also general enough to include the cases of two outgoing shock waves created 
by the focusing of compressive rarefaction waves (as simulated in |23]), or two 
incoming shock waves resulting in one outgoing shock wave and an outgoing 
rarefaction wave. 

Even though the Groah- Temple theory [7] establishes existence of C°'^ shock 
waves before and after interaction, and the work of Vogler numerically simulates 
the detailed structure of the metric at a point of shock wave interaction, the 
mathematical theory still lacks a complete proof that rigorously establishes the 
detailed structure of shock wave interactions summarized in Deflnition 13. II Such 
a proof would be very interesting, and remains to be done. 

4. Functions C°'^ Across a Hypersurface 

In this section we flrst deflne a function (or a metric) being C*^'^ across a 
hypersurface. We then study the relation of a metric being C°'^ across a hyper- 
surface to the Rankine Hugoniot jump condition through the Einstein equations 
and derive a set of three equations central to our methods in Sections |5] - [7) 
Finally we derive a canonical form functions C*^'^ across a hypersurface can be 
represented in, which is fundamental to the proof of our main theorem. 
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Definition 4.1. Let be a smooth hypersurface in some open set A/" C M"^ 
with a normal vector-field nowhere lightlike. We call a function f G C^'^{^J') 
"Lipschitz continuous across S ( or C"'"*^ across Ti), if f is smooth m A/'\ S and 
limits of derivatives of f to T, exist on each side of S separately and are smooth 
functions. We call a metric g^y Lipschitz continuous across S in coordinates 
if all metric components are C^'^ across S. 

For us, "smooth" means enough continuous derivatives so that regularity is 
not an issue. Usually / G C'^{M\ S) with left-/right-limits of first and second 
order derivatives existing suffices. 

The main point of the above definition is that we assume smoothness of / 
away and tangential to the hypersurface S, but allow for the normal derivative 
of / to be discontinuous, that is, 

UAn" 7^ 0, (4.1) 

where rf is normal to S with respect to some (Lorentz-) metric g^^^ defined on M. 
Moreover, the continuity of / across S implies the continuity of all derivatives 
of / tangent to S, i.e., 

[f,<rV = 0, (4.2) 

for all v'^ tangent to S. 



In the following we clarify the implications of Definition 4.1 , particularly (4.2 ) 



together with the Einstein equations on the RH conditions (3.4), (3.5). For 



this, consider a spherically symmetric spacetime metric in SSC (1.1) and assume 



the first three Einstein equations (2.21)-(2.23) hold and the stress tensor T is 



discontinuous across a smooth radial shock surface, described in the (t, r)-plane 
by 7(t) as in (3.1)-(3.3). To this end, condition (4.2) applied to each metric 



component g^i, in SSC, c.f. (2.20), reads 

m = 
[A] = 



-x[Br 



(4.3) 
(4.4) 



On the other hand, the first three Einstein equations in SSC (2.21 )-(2.23) imply 

[Br] = /tA5V[T°°], (4.5) 
[Bt] = -KAB'r[T^\ (4.6) 
[Ar] = KAB^r[T^^]. (4.7) 

Now 



using the jumps in Einstein equations (|4.5|)-(|4.7|), we find that (|4.3|) is 

is 



equivalent to the first RH condition (3.4), while the second condition (4.4 



ions 


4.5 


)-( 


4.7 


(2.21 


-( 


2.23 


). ' 



because At does not appear in the first 
tie result, then, is that in addition to the 
assumption that the metric be C°'^ across the shock surface in SSC, the RH 



conditions (3.4) and (3.5) together with the Einstein equations (4.5)-(4.7), yield 
only one additional condition]^ over and above (4.3) and (4.4), namely, 

[Ar] = -x[Bt] . (4.8) 



'^This observation is consistent with Lemma 9, page 286, of 
condition need to be imposed to meet the full RH relations. 



where only one jump 



REGULARITY SINGULARITIES 



17 



The RH conditions together with the Einstein equations will enter our method 
in Section Is] only through equations (4.8), (4.3) and (4.4). In particular, our 



proofs of Theorems |6. 1| and |7. 1| below will rely on the Einstein equation and the 
RH condition only in the sense of (4.8), (4.3) and (4.4). 

The following lemma provides a canonical form for any function / Lipschitz 
continuous across a single shock curve 7 in the [t, r)-plane, under the assumption 
that the vector n^, normal to 7, is obtained by raising the index in (3.3) with 
respect to a Lorentzian metric C°'^ across 7. 

Lemma 4.2. Suppose f is C^'^ across a smooth curve 7(t) = {t,x{t)) in the 
sense of Definition 4-1. t G (— e, e), in an open subs et M ofM?. Then there 



exists a function $ G (A/") such that 

1 



/(t,r) = -i^it) |x(t) -r| + <l>(t,r). 



if and only if 



G CH 



-e,e 



is a 1-form normal to v^{t) 



(4.9) 



(4.10) 



where n^{t) = (x(t),— 1) 
raised and lowered by a Lorentzian metric g^^^ which is C^'^ across 7 



'^^{t) and indices are 



Proof. We first prove the explicit expression for Lp in (4.10). Suppose there exists 

-- x{t) 



$ G C^i^M) satisfying ( [49| , defining X(t, r 
1 



with H{X) :-- 



- r this implies that 

-1 ifX<0 
+ 1 if X > ' 



(4.11) 



2 and X ^n'^ 



(4.12) 



where we use ^|X| = H{X) and [$,^] = 0. Since [H{X)] 
by (4.10), we conclude that 

Solving the above equation for ip we obtain the expression claimed in (4.10). The 
regularity of ip follows from our incoming assumption on the metric and /. 
We now prove the opposite direction. Suppose (p as defined in (4.10) is 
regular. To show the existence of $ G C^(A/') define 

1 



$ = |X| 



(4.13) 



then (4.9) holds and it remains to prove the regularity of $. It suffices to 
prove 



= 



since $ G C'^{Af\'-f) follows immediately from (4.13) and the regularity of / 
and (p away from 7. By assumption / satisfies (4.2), so that 



The right hand side vanishes since v'^ (t) = ^l,x(t)) and thus [<l>,^ 
Finally, ip defined in (4.10) together with 
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show that 

This completes the proof. □ 
In words, the canonical form (4.9) separates off the C^'^ kink of / across 7 into 



the function \x{t) — r\, from its more regular behavior incorporated into the 
functions ip, giving the strength of the jump, and $, which encodes the remaining 
behavior of /. 



In Section 7 we need a canonical form analogous to (4.9) for two shock curves 



which is provided by the next corollary. To this end, suppose timelike shock 



surfaces described in the (t, r)-plane by, ^{t) = {t,Xi{t)), such that (3.6) - (3.8) 
applies. To cover the generic case of shock wave interaction, we assume each 7j(t) 
is away from t = with the lower/upper-limit of the tangent vectors existing 



up to t = 0. For our main results (Theorem 7.1 and 1.1) it suffices to consider 
the upper {t > 0) or lower part {t < 0) of a shock wave interaction (at t = 0) 
separately, whichever part contains two shock waves that interact with distinct 
speeds. In the following we restrict without loss of generality to the lower part 
of a shock wave interactions, that is, to = {{t,r) : t < 0}. 

Corollary 4.3. Let •jiit) = (t, Xi(t)) be two smooth curves defined on I = (— e, 0), 
for some e > 0, such that the limits lim 7j(t) = (0,ro) and xAO) = lim xAt) 

both exist for i = 1,2. Let M be an open neighborhood of p = (0,ro) in and 
suppose f defined on A/'flM^ is in C°'^ across each curve 7j((— e, 0)). Then there 
exists a function ^ E (J\f fl ]R?_) , such that 

i=l,2 

for all (t, r) in M Pi if and only if 

<^,(t) = e C°'i(/), (4.15) 



where {ni)^{t) = {xi(t), —1) is a 1-form normal to v^it) = 'Jiit), (for i = 1,2), 
and indices are raised and lowered by a Lorentzian metric C^'^ across 7. 

Proof. The Corollary follows by the same arguments as in the proof of Lemma 



4.2 on each of the curves 7j, {i = 1,2), but defining 

(4.16) 



$ - , 

2 

i=l,2 



for Xi{t, r) := Xi{t) — r, in place of ( |4.13[ ) and using that H(Xi) is discontinuous 
across 7^, but [H{Xi)]i = for / 7^ i. (In particular, / meets ( 4.2[ ) across each of 



the curves 7^.) □ 



In Section [7| we use the canonical form (4.14) to characterize all Jacobians in 



the {t, r)-plane, that could possibly lift the metric regularity from C^'^ to C^'^ 
on a deleted neighborhood of the point p of shock wave interaction, unique up 
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to addition of a function $ G C^. It is precisely the regularity across 7 that 
forces the determinant of those Jacobians to vanish at p when taking the limit, 
resulting in a singular metric at p. 

5. A Necessary and Sufficient Condition for Smoothing Metrics 

In this section we derive a point-wise condition on the Jacobians of a coordi- 
nate transformation, necessary and sufficient for the Jacobian to lift the metric 
regularity from to C^'^ in a neighborhood of a point on a single shock surface 
S. This condition is the starting point for our strategy and lies at the heart of 
the proofs in Section |6] and [7j We begin with the covariant transformation law 

gap = J^Jp9t.u, (5.1) 

for the metric components at a point on a hypersurface S for a general C^'^ 
coordinate transformation x'^ — )■ x°', where, as customary, the indices indicate 
the coordinate system. Let Jj^ denote the Jacobian of the transformation, that 
is. 

Now, assume the metric components g^j^i, are only Lipschitz continuous across S, 



with respect to coordinates x^. Then, differentiating (|5.1|) with respect to ^ 



and taking the jump across S we obtain 

[9ap,,] = J^Jl[9,un] + 9,uJ'AJh] + 9,.Wl^^,] , (5.2) 
where [/] = fh — fn denotes the jump in the quantity / across the shock surface 



S, c.f. (2.15). Since both g^^ and are C°'^ across S, the jumps are only on 
the derivative-terms. Now, gap is in C^'^ if and only if 

[gapa] = (5.3) 



for every a, /3, 7 G {0, ...,3}, while (5.2) implies that (5.3) holds if and only if 



[■J^,]Jp9,u + [Jp,,]J!i9,u + €Jl[g,ua] = 0. (5.4) 

(5.4) is a necessary and sufficient condition for smoothing the metric from C°'^ 
to C^'^ by the means of a coordinate transformation, (c.f. Corollary 5.2). 



(5.2) also holds for Jacobians that are just Lipschitz continuous, possibly hav- 
ing discontinuous derivatives away from E, provided one defines the jump [■] 
along some (arbitrarily) prescribed curve, that is, define 

\u\ = lim u o c(t) — lim u o c(t) 

for some curve c transversal to S such that c(0) G S. 

Note, if the coordinate transformation is C^, the jumps in J^^ vanish, and 



(5.4) reduces to 

= JaJl3[gf_iu,'r], 

which is tensorial. Now, since the Jacobians are non-singular, we find that it is 



precisely the lack of covariance in (5.2) for C^'^ transformations, providing the 
necessary degrees of freedom, (namely [Ja^^]), making it possible for a Lipschitz 
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metric to be smoothed by coordinate transformation at points on a single shock 
surface. This illustrates, there is no hope of lifting the metric regularity within 

db C 3jt)lclS. 



Let's consider (5.4) as an inhomogeneous linear system with unknowns [J^ 



7J' 



interpreting the J^-factors as free parameters. A solution of (5.4) alone does not 
ensure the existence of a Jacobian, since ( 5.4[ ) does not yet ensure the existence 
of C°'^ functions that take on the values [J^^] and satisfy the integrability 



condition (2.4), necessary for integrating to coordinates. It is therefore crucial 
to impose in addition to (5.4) an appropriate integrability condition, namely 

[JU = [JU- (5.5) 



In the following we solve the linear system obtained from (5.4) and (5.5) for 



\Jaa\i subject to Lemma 



5.1 



To simplify (5.4), we subsequently restrict to 
spherically symmetric spacetimes and assume that gfj_u denotes the metric in 
Standard Schwarzschild coordinates]^ ( |2.20 ). To this end, suppose we are given 
a single radial shock surface S in SSC locally parameterized by 

E(t,^,0) = (t,x(t),^,0), (5.6) 

described in the (t, r)-plane by the corresponding shock curve 

7(t) = (t,x(t)). (5.7) 

For such a hypersurface in SSC, the angular variables play a passive role, and the 
essential issue regarding smoothing the metric components by C^'^ coordinate 
transformations, lies in the atlas of (t, r)-coordinate transformations. Thus we 
restrict to the atlas of (t, r)-coordinate transformations, which keep the SSC 



angular coordinates fixed, c.f. (2.20). Then 



( 4 


J\ 





\ 


Jo 


J{ 














1 













^1 



(5. 



Lemma 5.1 



with the coefficients just depending on the SSC t and r, which implies 

[J^^^] = whenever /i e {<^, oi a e {2,3} or /3 G {2,3}. (5.9) 

Now, ( |5.4[ ) and (5.5 ) form a linear inhomogeneous 8x8 system in eight unknowns 
The following lemma shows that this system is uniquely solvable. 

Consider a metric in SSC, 

g^^dx'^dx" = -A(t, r)df + B{t, r)dr^ + r^dQ^ , (5.10) 

let S denote a single radial hypersurface ( 5.6[ ) across which g^u is Lipschitz con- 
tinuous and assume (5.9). Then there exists a unique solution [J^o-] of the system 
(5.4) together with ( |5.5 ). The solution is given by: 



[Jo, 



[A 
A 



^0 



[A] 
A 



J'. 



[Ar] 

A 



jt , [JMjr 



Indices /z, v and a always refer to SSC, if nothing else is mentioned. 
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[Jit] 



A 
B 
B 





\AA 
A 


\ 




1 

~ ~2 


-4 + 


[Bt] 
B 


^^); 


[^0,r] 


1 

~ ~2 


-A + 


[Bt] 
B 




[J'lr] 


1 

~ ~2 



A 


-J{ + 


\BA 

A 


[Bt] 


4 + 


[Br] 


B 


B 


[Bt] 
B 


Jl + 


[Br] 

B 



J[ (5.11) 



{We use the notation /U, G {t,r} and a,/3 & {O?!}? that t,r are used to 
denote indices whenever they appear on the Jacobian J^.) 



Proof. We first setup the 8x8 linear system formed by (5.4) and (5.5) for tlie 

tlie (t, r)-plane, c.f. (5.9). 
are tliose for fi G {t,r} and a,f3 E {0,1}. Now, tlie 



metric (5.10) and the case of Jacobians in the (t, r)-plane, c.f. (5.9). The only 
nonzero terms of [J^^] 



integrability condition (5.5) reduces to 

[Jii] 

and [Jo^] 



[Jlo] 
[Jlo] 



(5.12) 



Using (5.12), we are left with six unknowns [J^^] and it remains to setup and 
solve (5.4) for the metric in SSC (5.10) and Jacobians in the (t,r)-plane. A 
major simplification of (5.9) is that the angular coefficients of the metric in SSC 

(5.13) 



(5.10) drop out and (5.4) reduces to the following 6x6 linear system 

Av 



with unknowns 



w 



[Jli]) 



where T denotes the transpose. The right hand side in (5.13) is given by 



for 



W = [Wi,W2,W3,W4,W5,W6) 

= '^{wooo, wolo, Who, wqoi, wqu, wm) 



(5.14) 

(5.15) 
(5.16) 



The matrix A reads 



A 



( 2a 
c 







V 





a 
2c 
2a 
c 









a 
2c 



2d 
b 








d 
2b 
2d 
b 




\ 






d 

2b ) 



(5.17) 



with a 



—JqA, 



b := JlB, c := —J^A and d := JqB. Observe the we already 
incorporated the integrability condition (5.12) into the linear system (5.13) and 
thus only considered 6 unknowns in (5.14). 
(5.13) has a unique solution given by, 

. , 2wod — bw^ — dwA 
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[Jli] 

[Jll] 



dw3 — bW4 

2AB\J\ '■ 
dwQ + 



2AB\J\ 
aw4 — W2a + cwi 

2AB\J\ ' 
cWi — aw3 

2AB\J\ ' 
2cu'5 — awQ 



CW3 



2AB\J\ 



(5.18) 



where \J\ = det (J^) = JqJ[ — J\Jq denotes the determinant of the Jacobian. 
It remains to change coordinates to SSC on the left hand side of (5.18 ), in order 



to recover the expressions in ( 5.11[ )j | From the chain rule and the continuity of 
the inverse Jacobian J°" we get 



and (for cr = t, r) 



I r 



Using the definition of Wi in (5.16) together with (5.18) and the jumps in the 



integrability condition (5.5), a straightforward (but lengthy) calculation leads to 

dslTl). 



In order to exemplify the procedure we outline the calculation leading to the 
expression for [Jq J in (5.11): 



[Jo,t 



[^0,0^ 



\J\ 



Ji 



(5.19) 



where we used Cramer's rule in the last equality. Now, substituting (5.18) to- 



gether with the definition of Wi in (5.16) and the integrability condition (5.5) a 
straightforward calculation leads to 



A 



■Jo + 



[A] 
A 



V 



as claimed in (5.11). Performing similar computations for the remaining expres- 
sion finishes the proof. □ 



By Lemma 5.1, (5.11) is a necessary and sufficient condition for [g, 



since it solves (5.4) and (5.5) uniquely. Putting (5.11) in words, the Jacobian 



0, 



must mirror the metric regularity, compensating for all discontinuous first order 
metric derivatives, to ensure [(70/3,7] = In fact, assuming the [J^,^]-terms 

" " ^" 'dsTTl) is 



to come from an actual Jacobian of a coordinate transformation, J^, 



^We want to consider the shock curve 7 in the coordinates we start in (SSC) and thus 
express the condition (5.11) in terms of the given coordinates x^^ . 



"'^"This is why we expect a C^'^ atlas to be generic for possibly lifting the metric regularity 
from C"'i to C^'i. 
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necessary and sufficient for raising the metric regularity to C ' at a point on a 



single shock surface, as is proven in the next Corollary 5.2 



Corollary 5.2. Let p be a point on a single smooth shock c urve 7, and let g^^ 

he the metric in SSC, C°'^ across 7 in the sense of Definition 4-i Let be the 



Jacobian of an actual coordinate transformation defined on a neighborhood M of 
p and assume is C^'^ across 7. Then the metric in the new coordinates, gap, 



is in C^'^{N) if and only if satisfies (5.11). 



Suppose there exist 

such that the metric in the new coordinates gas is C^'^ regular, 



Proof. We first prove that ga 
coordinates x' 
then 



e Ci'i(Ar) implies (5.4) 



V a, /3, 7 e {0,...,3}. 



This directly implies (5.4) and since are the Jacobians of an actual coordinate 



5.1 



transformation they satisfy the integrability condition (5.5) as well. By Lemma 
the jumps in the derivatives of the Jacobian [J^,^,] then satisfy (5.11). 



We now prove the opposite direction. Suppose the Jacobians satisfy (5.11) 



by Lemma 5.1, they then meet the smoothing condition (5.4), which implies all 
metric derivatives gaj3,'y to match continuously across the shock curve 7, that 
is, [(70/3,7] = for all a, /3, 7 e {0, ...,3}. Since g^j_u and are assumed to be 
smooth away from 7, it follows that gap € C^{M). Moreover, the existence of 
the limit towards 7 of all ffist and second order metric derivatives is postulated, 
(c.f. Definition 4.1), which then implies gap G C^'^{N'). Namely, for any g G S 
one can bound 

\9cp,^{t, r) - gapM \<c\{t,r)- (g°, g^) | , (5.20) 
by a finite constant c, for instance choose 

c= V ( sup |(^«/3,^<,)^(g)| + sup |(^«/3,^<,)^(g')| )+ y] sup |^„;3,^(t, r)| , 

a,p,'y,S ^1'^^' ^ a,p,-y (^'-)^^' 

=0,1,2,3 =0,1,2,3 

for A/"' C A/" a compact neighborhood of q and S' = S fl A/"'. □ 



Lemma 5J^ and Corollary 5^ establish the remarkable result that there is no 
algebraic obstruction to lifting the metric regularity. Namely, the smoothing 



condition (5.4) has the solution (5.11), and its solvability is neither connected to 
the RH conditions (3.4)-(3.5) nor the Einstein equations. However, by flEl [TO] , 



the RH conditions must be imposed for coordinate transformations smoothing 
out the metric to exist. The point here is, that to prove the existence of coordi- 
nate transformations lifting the regularity of SSC metrics to C^'^ at p G S, one 
must prove that there exists a set of functions Jj^ defined in a neighborhood of 
p, such that (5.11) holds at p, and such that the integrability condition (2.4), 



necessary for Jj^ to be the Jacobian of a coordinate transformation, holds in a 
neighborhood of p. In Section [6} we give an alternative proof of Israel's Theorem 
by showing that such exist in a neighborhood of a point p on a single shock 
surface if and only if the Einstein equations and the RH condition hold. Such 
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a Jacobian fails to exist if two radial shock surfaces intersect, as we prove in 
Section [71 

6. Metric Smoothing on Single Shock Surfaces and a 
Constructive Proof of Israel's Theorem in Spherical Symmetry 



We have shown in Corollary 5.2 that (5.11) is a necessary and sufficient con- 
dition on [Jj^f^] for lifting the SSC metric regularity to C^'^ in a neighborhood of 
a shock curve, provided the value [Jao] comes from an actual Jacobian J^. We 
now address the issue of how to obtain such Jacobians, that is, we study how 



to construct a set of functions satisfying (5.11) on the shock curve, while the 
integrability condition (2.4) holds in a neighborhood of the curve. For a single 



shock surfaces, this can be done if and only if the RH condition hold, as recorded 
in the following theorem. 

Theorem 6.1. (Israel's Theorem) Suppose g^j^u is an SSC metric that is C^'^ 



across a radial shock surface S in the sense of Definition ^.i, such that it solves 



the Einstein equations (2.21) - (2.24) strongly away from E. Assume T^'^ is 
smooth away from S and let p be a point on S. Then the following is equivalent: 
(i) There exists a coordinate transformation of the {t,r) -plane, defined in a 
neighborhood JV of p and with its Jacobian J^^ being C^'^ across S, such 
that the transformed metric components gap are C^'^ functions of the new 
coordinates. 



(a) The RH conditions (3.4) - (3.5) hold on in a neighborhood of p 



Theorem is the spherically symmetric version of results first proven by 
Israel [10], see also [18]. In the remainder of this section we develop a new 



method to prove Theorem 6.1 , which also allows us to address shock interactions 
in Section [7| finally leading to our main result. Theorem 1.1 The first step 



is to construct Jacobians acting on the (t, r)-plane that satisfy the smoothing 
condition (5.11) on the shock curve, the condition guaranteeing [5'a/3,7] = 0. For 



a single shock curves, both the RH conditions and the Einstein equations (in the 



form of (4.8)) are necessary and sufficient for such functions Jj^ to exist, c.f. the 



discussion in the end of Section [5} as presented in the following lemma. 
Lemma 6.2. Let p be a point on a single shock curve 7 across which the SSC 



metric g^j^i, is Lipschitz continuous in the sense of Definition 4-1 Then there 
exists a set of functions defined on M being C°'^ across S and satisfying 
the smoothing condition (5.11) on ^ r\ M if and only if (4.8) holds on 7 fl M . 
Furthermore, all that are C^'^ across S and satisfy (5.11) on'-fClAf are given 
by 



4{t,r) 
Jiit,r) 



AA o 7(t) 
4:A o 7(t) 

[Bt]<f){t) + [BrMt) 

AB o 7(t) 



\x{t) -r\ + r) 
x{t) -r\ +N{t,r) 
\x{t) - r| + Q{t,r) 
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J[(t,r) 



[B,]u{t) + [BrUt) 



AB o 7(t) 

for arbitrary functions $, fi, Z, G C^^M) and 



— r| + r) 



Z O J . 



(6.1) 



(6.2) 



Proof. Suppose there exists a set of functions satisfying (5.11) and being C°'^ 
across E, which implies 

Kt] = -i[JL] (6.3) 



for all /X G {t, t} and a G {0, 1}. Combining (6.3) for the special case /i = t and 
a = with the right hand side in ( |5.11 ) leads to 



A 



'0 



A 



A 



Using now the jump relations for the metric tensor, (4.3) - (4.4), finally gives 
[Ar] = -x[Bt], that is, (|48|. 



For proving the opposite direction it suffices to show that all t and r derivatives 



of J^, defined in the above ansatz (6.1), satisfy (5.11) for all /x G {t,r} and 
a G {0, 1}. Since (6.2) implies the identities 

= 407, z/ = J{ o 7, a; = Jo" o 7, ( = J[ o 7 , (6.4) 

and since the regularity of $, Q, Z and implies 

[f/^] = = [f/i], for f/ = $, Z or AT, (6.5) 

(4.3), (4.4) and ( |4.8 ) immediately show that the Jacobian ansatz (6.1), satisfies 
(5.11). This proves the existence of functions satisfying (5.11). 



For the supplement, applying Corollary |4.2| confirms that all functions satis- 

□ 



fying (5.11) can be written in the canonical form (6.1) 



To complete the proof of Israel's Theorem, it remains to study under what 



condition the functions defined in (6.1 ) satisfy the integrability condition (2.4) 
in a neighborhood of the shock, ensuring that is integrable to coordinates. 
This is accomplished in the following two lemmas. The first lemma gives an 



equivalent form of the integrability condition (2.4), adapted to the free functions 



$, f2, Z and of the Jacobian ansatz (6.1). The main step in its proof is to 



write (2.4) in SSC. 



Lemma 6.3. There exist coordinates x" such that Jj^ defined in (6.1) satisfy the 



integrability condition (2.4) if and only if the functions ^,il,N and Z, in (6.1), 
satisfy the following system of PDE's: 



{a\X\ + $4) (/3|X| +N) + ^r {e\X\ + Z) - {a\X\ + $) (^/3\X\ + Nt) (6.6) 

-Nr{6\X\+n) + fH{X) = 



6\X\+nt) {f3\X\+N) + nr (e|X| 



Z) - (e|X| + Zt)(a|X| + $) (6.7) 
-Zr i6\X\ + hH{X) = , 
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where X(t, r) 



x(t) — r, H{X) denotes the Heaviside step function, 

[Ar](j){t) + [Bt]uj{t)_ 



a{t) 

m 
m 

e{t) 



AA o 7(t) 

[AMt) + [Btm _ 

4A o 7(t) 

[Bt](P{t) + [BMt) . 
AB o 7(t) 

[BMt) + [Br]m _ 

AB o 7(t) 



(6.8) 



and 



f = {(35- ae) \X\ + axN - + (5VL - aZ 
h = {(36 - ae)x\X\+ 6xN - ex^ + eQ - 6Z . 



(6.9) 



Proof. Suppose there exist coordinates x", such that satisfies the integrabihty 



condition (2.4), that is 



From the chain rule we find that (6.10) imphes 



(6.10) 



(6.11) 



(where x^ denote SSC). Substituting the Jacobian defined in (6.1) into (6.11) 



proves that $,n,A^ and Z satisfy (6.6)-(6.7). 



We now prove the opposite direction. Observe that, in hght of (6.1), the 



integrabihty condition (6.11) is in fact equivalent to (6.6)-(6.7). It remains to 



prove that (6.11) implies the existence of coordinates x° such that (6.10) holds. 



For this, we show that (6.11) implies the integrabihty condition on the inverse 



Jacobian J^, that is. 



ja 



(6.12) 



For the moment, assume (6.12) holds, then the function are integrable to 
coordinates x", (c.f. Lemma A.l), and the resulting function (x^) i— )■ (x") is 



bijective, which allows us to introduce 



Jit 



dx^ 
dx°' 



(6.13) 



By the chain rule, (6.13) gives the point-wise (linear algebraic) inverse of J° 



Moreover, from (6.13), the satisfy the integrabihty condition 6.10 by the 



once we have shown that (6.11) imply (6.12). 



We now prove this implication. Suppose 6.11 holds, which is equivalent to 



commutativity of partial derivatives. We conclude that Lemma 6.3 is proven 



jr jt _ jr jt 
jr jr jr jr 



jt jt _ jt J. 

p r _ /* / 



t 

r 

0,t 



(6.14) 

(6.15) 
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A direct computation, using (6.14) - (6.15), verifies (6.12). In more detail, the 



linear algebraic inverse of is given by 

J? 
J} 



't 



Jl 



1 

\J\ 



jr _ Jt 



- P' 
^0 



7* 

^0 



(6.16) 



where \J\ denotes the determinant of J^. From (6.16), a straightforward com- 
putation gives 



7° 



jr 

( f f _ f r ] _^ ( r p _ / p ) 



Exchanging the first term using (6.15) and the second term using ( |6.14 ), we get 



A 

|J|2 


r,t 1 



( p r - 7* r ) - iP P - 7* 7* 



J 



where (6.16) enters the last equality. A similar computation verifies the remain- 



ing equation in (6.12), 



t,r 



n 



This shows that (6.11) implies (6.12) and 

□ 



completes the proof. 

The proof of Israel's Theorem is complete once we prove the existence of 
regular solutions ^,Q,N and Z of (6.6) and (6.7). For this, it suffices to choose 
N and Z arbitrarily, then (6.6) - (6.7) reduce to a system of 2 linear first order 
PDE's for the unknown functions $ and il. Thus ( |6.6 ) - (6.7) can be solved along 
characteristic lines and the only obstacle to solutions $ and Q with the necessary 
regularity is the presence of the (discontinuous) Heaviside function H{X) in 
(6.6) - (6.7). However, the coefficients of H{X), f and h, vanish precisely when 
the RH jump conditions hold on 7, as stated in the next lemma. 

Lemma 6.4. Assume the SSC metric g^i, is across 7, in the sense of Def- 



inition and solves the first three Einstein equations strongly away from 7. 

(6.17) 



Let f and h be the functions defined in (6.9). Then, 



= h 



7 



if and only if the RH conditions (3.4) - (3.5) hold on 7. 



Proof. The first terms of / o 7 and /i o 7 in (6.9 ) drop out, since X o 7 = 0. Now, 
using (6.2) gives 



/07 
hey 



axu — f3x(f) + (3oj — a( 
5xu — ex(f) + ecu — 6(. 



(6.18) 



In (|6.18|), replace a and /3 by their definition, (6.8), then a straightforward com- 

/o7 = 

i[Ar]+x[Bt]){(l)C-uu) = 0. (6.19) 



putation shows that 
is equivalent to 
Now, using 



(0C-z/u;)=det(J^ 07)^0, 



(6.20) 
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we conclude that / o 7 = if and only if (4.8) holds, which is equivalent to the 



second RH condition (3.5). 



Similarly, replacing 5 and e in (6.18) by their definition, (6.8), a straightforward 







computation shows the equivalence of 

/l o 7 
and 

{[Bt]+x[Br\){^C-uuj)=^. (6.21) 

Now, using again ( [6^ , it follows that /i o 7 = if and only if (Q holds 
which is equivalent to the first RH condition (3.4). This completes the proof. 



□ 



In summary. Lemma 6.4 asserts that $, f2 G C^{Af) if and only if the Rankine 



Hugoniot conditions (3.4)-(3.5) hold. We are now in the position to prove Israel's 
Theorem: 



Proof, (of Theorem 6.1) Assume there exist coordinates x'^ 



We now prove that the RH condition (3.4)-(3.5) hold as a consequence. 



such that Qap is in 



Recall that the first RH condition is implied by the continuity of the metric 
see the discussion following (4.3). It remains to verify (4.8). The C^'^ regularity 
of Qafi implies that the Jacobian satisfies (5.11) on 7. Thus, by Lemma 6.2 



the Jacobians are of the canonical form (6.1) with $, f2, and Z being C 



regular, hence satisfying (6.5). Moreover, the functions $, fi, and Z satisfy 
the integrability condition. Taking the jumps of the integrability condition (6.6)- 

and /i o 7 



(6.7) and using (6.5) implies /07 



0. By Lemma 6.4[ we conclude 
This proves 



that (4.8) holds and therefore also the second RH condition (3.5) 

implies 



following from (i) 
We now show that 



Suppose the RH conditions (3.4)-(3.5) and 



thus (4.8) hold on 7 fl A/". By Lemma 6.2, the functions defined in (6.1) 



satisfy the smoothing condition (5.11) and thus lift the metric regularity from 
C°'^ to C^'^. By Lemma 6.3, to ensure the existence of coordinates such that 
Jit 



1^, we need to prove the existence of regular functions $, Vt, N and Z, 



that satisfy the integrability condition (6.6). For this, without loss of generality 



we choose A^ and Z arbitrarily and solve the resulting linear system of PDEs, 
(Q, for $ and 



In more detail, choose A^ and Z to be smooth functions on M subject to the 
condition 

C^xu (6.22) 

on the shock curve, where z/ = A^ o 7 ^ and C = Z o 7 ^ 0. Due to ( |6.22[ ), the 
characteristic curves of (6.6) are transversal to the shock curv^^ ^{t) = {t, x{t)). 
This allows us to prescribe initial data on the shock curve, i.e., to prescribe 
$07 = and n o 7 = cj, then the derivative-terms in (6.6), a and 6, are 



given functions. Now, (6.6) is a linear symmetric hyperbolic system of first order 



^""^Due to the metric regularity, (4.4) always holds, c.f. Section |4j 

^^This holds, since {v, C) is tangent to the characteristic curve on 7, c.f. fTT] for an intro- 
duction to the method of characteristics. 
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PDEs for $ and ^l. Choosing smooth initial data, the method of characteristics 
in [TT], (chapter 2.5; exchanging (t,x) by (r, t)), shows existence of a solution 



($, f2) on some open neighborhood of 7 as follows: Integrating (6.6) along its 
characteristic curves, the method in [TT] yields a solution to the resulting 
integral equation, (even for the discontinuous coefficients of lower order terms in 



(6.6)-(6.7)). Once we prove the regularity of this solution it is clear that it 



also solves the original PDE (6.6). 



We now prove the regularity of the solution ($, Q). Away from the shock 



curve 7 (and tangential to 7) all coefficients in (6.6) are smooth, thus we can 



follow the method in jTT] proving smoothness of ($,fi) away from the shock 
curve. In particular, choosing initial data, and Z such that the limit of first 
and second derivatives exists on each side of 7, then $ and fl are C°'^ across 7. It 



remains to prove a regularity across the shock curve, that is, (6.5). For this 



the vanishing of / and h in (6.17), the continuity of ($, f2) and the integrability 
condition (6.6) imply that derivatives of ($, i^) along the characteristic curves 
are continuous. That is, 



[$i]z/+[$,]C 




0. 



(6.23) 



(Note, (i^, C) is tangent to the characteristic curve of (6.6) at 7.) On the other 
hand, ($,fi) being C°'^ across 7 implies (c.f. ( 4.2[ )) 

[ilt] + [ar]x = 0. (6.24) 



Now, since (6.22) ensures that the shock curve is non-characteristic, (6.23) and 



(6.24) are independent condition, which then yield the desired regularity 



across 7, (6.5). We conclude that ($,^^) are regular on some neighborhood 
A/" and smooth on A/" \ 7. Therefore, is across the shock curve. 

Moreover, we can choose initial data and the free functions, 0, C, such that 
det (J) 7^ everywhere on the shock curve and by continuity of the solution 
det (J) is non- vanishing in some neighborhood of 7. 

In summary, assuming the RH conditions on 7, the Jacobian is C^'^ across 



the shock, satisfies the smoothing condition (5.11) on 7 and satisfies the integra- 



bility condition on some neighborhood JV. Therefore, by Corollary 5.2, maps 
the SSC metric, g^u, to a C^'^ regular metric gajs- This completes the proof of 
Theorem [6l □ 

In fact, in Theorem|6| one could replace the assumption of being C°'^ across 
E by the assumption that the coordinate transformation is C^'^ regular. Namely, 
if there exist a C^'^ regular coordinate transformation smoothing the SSC metric 
to C^'^, one can assume without loss of generality that the Jacobians are C°'^ 
across S, since the smoothing condition (5.1) is a point-wise condition on the 
shock surface. Note, (5.1) is also valid for Lipschitz continuous Jacobians not 
smooth away from the shock surface, provided one defines the jump [■] in terms 
of some curve transversal to S, c.f. Section [s] Moreover, proving that (i) implies 
(ii) is not affected, since the Jacobian constructed is C°'^ across 7 anyways. 
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7. Shock Wave Interactions are Regularity Singularities in GR; 
Transformations in the (t, r)-PLANE 



The main step in the proof of Theorem l.l| is to show the result for the smaller 
atlas of coordinate transformations in the {t, r)-plane first. That is, to prove there 
do not exist C^'^ coordinate transformations of the (t, r)-plane (i.e., keeping the 
SSC angular variables fixed) lifting the regularity of the SSC metric g^i, from 
C^'^ to C^'^ in a neighborhood of a point p of regular shock wave interaction in 
SSC. We then prove in Section |8] that no such transformation can exist within 
the full C^'^ atlas that transforms all four variables of spacetime, including the 



SSC angular variables, leading to Theorem 1.1 We formulate the main step for 
lower shock wave interactions in in the following theorem, which is the topic 
of this section. A corresponding result applies to upper shock wave interactions 
in Ml, as well as to shock wave interactions in a neighborhood of p. 



Theorem 7.1. Suppose that p is a point of regular shock wave interaction in 
SSC for the SSC metric g^y, in the sense of Definition 3.1. Then there does not 
exist a C^'^ coordinate transformation o [x^)^^ of the {t,r)-plane, defined on 
M n Ml for a neighborhood M of p in M^, such that the metric components gap 
are functions of the coordinates x"" in Af (1 Ml and such that the metric has 
a non-vanishing determinant at p, {that is, such that limdet {gapiQ)) 7^ 0). 



Theorem lA implies the non-existence of coordinates on a neighborhood TV of 
j9 in M^. In the remainder of this section we present the proof of Theorem 7.1 



which mirrors the constructive proof of Theorem 6.1 in that it uses the extension 



of our Jacobian ansatz (6.1 ) to the case of two interacting shock waves. But now. 



to prove non-existence, we must show the ansatz is general enough to include all 
C^'^ Jacobians that could possibly lift the regularity of the metric. We conclude 
the proof by showing that this canonical form is inconsistent with the assumption 
that det (ga/s) 7^ at p, by using the continuity of the Jacobians up to p. 

To implement these ideas, the main step is to show that the canonical form 



(4.14) of Corollary 4.3 can be applied to the Jacobians in the presence of a 



shock wave interaction. The result is recorded in the following lemma: 

Lemma 7.2. Let p be a point of regular shock wave interaction in SSC in the 

corresponding to the SSC metric g^y defined on M r\ Ml . 



3.1 



sense of Definition 
Then the following is equivalent: 

(i) There exists a set of functions J^, defined on Af nM?_, which is C 



0,1 



across 



7j n A/" and satisfies the smoothing condition (5.11 ) on 7^ fl A/", for i = 1, 2. 

1,2. 



(a) The RH condition (4.8) holds on each shock curve 7^ fl Af, for i 
Moreover, if (i) or (ii) is valid, then the functions assume the canonical form 

Joit, r) = Yl I^^W - r| + $(t, r), 

i 
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•^0 r) = Yl W - '"I + '^)' 

i 

Jlit, r) = Y^ \xi{t) -r\ + Z{t, r) , 



(7.1) 



where ^,n,Z,N e C^AT n 

ai(t) 

m 



and 






(j)i{t) + [Bt], ujiit) 




4Ao7,(t) 




Ui{t) + [Bt]i Ci{t) 




4Ao7,(t) 




(j)i{t) + [Br]i Ui{t) 




4So7,,(t) 




Ui{t) + [Br]i Ci{t) 


45o7i(t) 



with 



$o7i, uJi = no'ji, = Zo7i, u^ = No-ii 



(7.2) 
(7.3) 



Proof. The proof is analogous to the proof of the single shock version in Lemma 
6.2 In more detail, suppose that there exist functions J/^ that are C'''^ across 
each 7i and meet the smoothing condition (5.11) on each 7,, for i 
being across 7^ implies (c.f. (4.2) and (6.3)) 



1,2. 



4^ 



[Kt]^ = -iAJL]^ for z = l,2. 



(7.4) 



Substituting into (7.4) the expressions for [JaJ\i from the smoothing condition 



(| 5.11 ) (with respect to 7i) and using that g^y is C°'^ across 7i, in the sense of 
|4.3[)-(|44)), finally yields ( [48| on 7^. Now, ( [X3| together with ( [48| imply the 
RH conditions on 7^ for z = 1,2. 

We now prove the implication from (ii) to (i). Suppose that the RH conditions 
hold on each 7j, i = 1,2, then (4.3), (4.4) and (4.8) hold on each shock curve. 



due to g^y being C°'^ across each 7j. Now, a straightforward computation using 
(4.3), (4.4) and (4.8) to eliminate Xi shows that the Jacobian ansatz defined in 



the existence of functions satisfying (5.11). 



(7.1 ) satisfies the smoothing condition (5.11 ) on each shock curve 7^. This proves 



From Corollary 4.3 we conclude that if functions as in (i) exist, then there 
exist functions $, f2, Z, iV G C^{Af), such that the can be written as in 
(7.1). This completes the proof. □ 



We now have a canonical form for all functions that meet the necessary 
and sufficient condition for [5^0/3,7] = 0, (5.11). However, for to be proper 



Jacobians, integrable to a coordinate system, we must use the free functions 



$, fl, Z and to meet the integrability condition (2.4). In the spirit of Lemma 



6.3, the following lemma gives a form of the integrability condition tailored to 



our Jacobian ansatz (7.1). To prove Theorem 7.1, the precise form of the in- 



tegrability condition is not relevant, except the structure of the coefficients of 
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the (discontinuous) H{X)-teTms, f and h. The main difference to (6.6) is the 
appearance of additional mixed terms in / and h, recorded in (7.6). 



Lemma 7.3. For the ansatz (7.1) the integrability condition (2.4) is of the form 
G(ipt, ipr, '-p, t, r) +^ FiH{Xi) = 0, 



(7.5) 



=1,2 



where Xi{t,r) = Xi{t) — r and H{.) denotes the Heaviside function, i.e., H{X) = 
-1 for X < and H{X) = 1 for X > 0. <^ ($, Q, Z, N) and ipi = 
^(<^o7j)(t). G is a (2-D) vector valued function continuously depending on 
its arguments and Fi ='^{fi, hi), where fi and hi are real valued functions con- 
tinuously depending on t and r, with their restriction to the shock curves given 
by 



fi°li 



hi 



li 



{oLiXiVi - l3iXi(j)i + iSiUi - aid) 
+ {aiXif3i - l3iXiai + - a^e/) | 
{5iXiVi - eiXi(f)i + eiUJi - did) 
+ {5iXi(3i - eiXiUi + ei5i - 6iei) \x. 



- xi( 



(7.6) 



for i = 1,2 and I ^ i. 



Proof. The proof is similar to the one of Lemma 6.3| in that the equivalence of 
the integrability condition (2.4) and the integrability condition in SSC (6.11), 
that is, 

follows by exactly the same arguments. (In fact, this equivalence holds inde- 



pendently of the Jacobian ansatz (6.1)). Substituting now the Jacobian ansatz 



for the intersecting shock curves (7.1) into (6.11) and separating continuous and 
discontinuous terms (i.e. terms containing H{X)) leads to (7.5). □ 



For proving Theorem 7.1, the point of interest in (7.5) lies in the additional 



mixed terms appearing in the coefficients fi and hi of the if(X)-terms. In con- 



trast to the single shock case (6.6), those mixed terms do not vanish on the shock 



curves 7^ by the RH condition alone. But, as a consequence of the regularity 
of $, Vt, Z and A^, fi and hi must vanish on the shock curves, which imposes an 
additional constraint. Taking the limit of this constraint to the point p of shock 
wave interaction finally yields the condition that det ( J^) vanishes at p, which 
immediately implies det {gajs) \ =0. 

To 



The main step in the proof of Theorem 7.1 is recorded in Lemma 7.4 



appreciate its significance, let's first compute the determinant of the Jacobian at 
the point of shock interaction. From ( |7.3[ ) together with ( 7.1[ ), we find 

det ( o 7,(t)) = ( J* J[ - J^) = mm - u,{t)u,{t) (7.7) 

and taking the limit to the point of shock wave interaction gives 



limdet (J^o74t)) 



0oCo - ^0(^0, 



(7. 
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where we use the notation 



00 = 



uq = hmcji(t) , Co = hm6(t) and z/q = hmz/j(t) . (7.9) 

t/'O t/'O t/'O 



Note, that the hmit in (7.9) is independent of the shock curve 7^ 
is assumed to be continuous and since the hmits 



1,2), since 



are assumed to exist, c.f. Definition 
lemma for the proof of Theorem |7.1 



3.1 From (7.8), the importance of the next 
is obvious. 



Lemma 7.4. Let p E M he a point of regular shock wave interaction in SSC in 
the sense of Definition 3.1 Assume the functions defined in (7.1) satisfy the 
integrability condition 

(7.10) 



on M r\ then the following equation holds at t = 0, 
1 



AB 



X1X2 1 



[Br]l[Br]2 {Xl - X2) (0oCo " J^O^^o) 



0. (7.11) 

(The metric coefficients A and B in (7.11) are evaluated at p = (0,ro) and 
(f)Q,...,Uo are defined in (7.9).) 



Proof. Suppose the functions J^, defined in our Jacobian ansatz (7.1), solve the 
integrability condition (7.10). Then, by Lemma 7.3, the functions $, Q, Z and 
satisfy the integrability condition in SSC (7.5). Taking the jump of (7.5) across 
7i gives 

G {ipt, ipr, ^, ipu t, r)] + V [F, H{Xj)]^ = 0. (7.12) 

i=i,2 

Using $, Q, Z,N G C^{M fl M.'^) and using that G and Fj depend continuously 
on their arguments, ( |7.12[ ) reduces to 



F, [H{X,)]^ = 0. 



Using further that [H{Xj)]i = 26ij, we get 

/o7,(t) = = /io7,(t), 

for all t G (— e, 0) and for i = 1,2. 

In the following we compute the limit of the equation /i o 7j(t) = as t 
approaches oj^ which directly leads to ( |7.11 ) and completes the proof. Without 
loss of generality we choose i = 1 (and thus / = 2 in (7.6)). Using (4.3) with 
respect to 71, which holds since the metric is across 71, the first term in 



"'^■'Note, one can prove that ho^i{t) = implies fo^i{t) ~ 0. The inverse imphcation holds, 
provided ±i{t) ^ 0. However, computing the limit of / o 7j(t) — a.s t approaches leads to 
( |7.11[ ) multiplied with ii(0), which does not suffice to prove Theorem |7.f [ unless we assume 
i.(0) ^0. 



34 



M. REINTJES 



(7.6) vanishes on 71 by the same arguments as in the proof of Lemma 6.4 We 
obtain 

/ii o 7i = - eiii«2 + £1^2 - ^162) \xi{.) - a;2(.)| . (7.13) 

For t 7^ 0, the initial assumption Xi(0) 7^ ^2(0) implies \xi{t) — X2{t)\ 7^ for all 
t sufficiently close to 0. Thus ho '-fi(t) = is equivalent to 

6iXil32 - exXxOL2 + ^1^2 - 5ie2 = 0, (7.14) 



sufficiently close to t = 0. By continuity, (7.14) holds at t = as well 



We next compute the limit of (7.14) as t approaches from below. From the 



definition of ctj and 5j in (7.2), together with the RH condition in the form (4.3) 



and (4.8), we get the identity 



•^i ^ '^i ; 



(7.15) 



for i = 1,2, where A and B are evaluated at p = (0, tq). Now, from the definition 
of €{ and Pi in (7.2), as well as (7.15) and the continuity of the Jacobian at p, we 
conclude that ( 7.14[ ) at t = is equivalent to 



+ 



. [Arh 


[Bth] 




\ [Ar], 


\Bt\i 


A 


B J 




I A 


B 


■ [Bth 
A 


[Brh 1 




\ 

A 


\Br\l 


B J 


B 



Co = 0. (7.16) 



Using (4.3) and (4.8) to eliminate [Bt]i and [Ar]i we find that the coefficients to 
61 and 62 in (7.16) are related as follows: 



Xl 
Xl 



[BtU_[Br 

A 

[Ar] 



B 

[Bt]r 



X1X2 



X^ 



1 

1 

V A +5 



A 
f X1X2 



[Br 
[Br 



A B 

for i = 1,2 and / 7^ i. Substituting (7.17) back into ( 7.16[ ) gives 

X1X2 !_ 
B 



(7.17) 



A 



[MBrhSi - x.lBrW Uo+{-[Br]26l + ^162) (o) = 0. (7.18 



The definition of 6i in (7.2) gives, at t = 0, 

[Bt]^ 00 + [Br 



AB 



for i = 1,2. Substituting the above equation into (7.18) proves that (7.14) is 
equivalent to (7.11). This completes the proof. □ 



Proof, (of Theorem 7.1 ) Assume there exists a C^'^-coordinate transformation in 
the {t, r)-plane mapping SSC to coordinates such that gap is in C^(A/'n]R?_) for 
some neighborhood M of p. Thus the Jacobian of the coordinate transformation 



satisfies the smoothing condition (5.4). Since (5.4) is a point-wise condition on 



the shock cures, also valid for Lipschitz continuous Jacobians not smooth away 
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from the shocks, we assume without loss of generahty that the Jacobians are 
C^'^ across each ji, {i = 1,2). 

Now, by Lemma 5.1, the Jacobian of the coordinate transformation sat- 
isfies the smoothing condition (5.11) on each of the shock curves ji, [i = 1,2). 

the Jacobians JH can be written in terms of 



Consequently, applying Lemma 



7.2 



the canonical form (7.1) with regular functions and Z. 



By assumption, the are integrable to coordinates and thus satisfy the in- 
tegrability condition (2.4), (c.f. Appendix lAl) . This together with being of 



the canonical form (7.1), allows us to apply Lemma 7.4, which yields that (7.11) 



holds at the point p of shock wave interaction, that is. 



1 / X1X2 



1 

~B 



[B,]l[Br]2 (Xl - X2) ( 0oCo - l^O^O , 



0. 



(7.19) 



To finish the proof of Theorem 7.1, observe that the first three terms in (7.19) 



are nonzero by our assumption that shock curves are non-null, and have distinct 



speeds at t = 0. In more detail, by our assumptions in definition 3.1 

[xi - ±2) ^ 0, 

while 

XiXi 1 , 

for time-like shock curves, since such curves satisfy Bx] < A. Finally, 



since otherwise [T 



[Br]i ^ 0, 

for all /i, z/ = 0,1, due to the Einstein equations 



(4.5)-(4.7), in contradiction to our assumptions in Definition 3.1 We conclude 

0oCo - ^^0^0 = (7.20) 
and using the explicit expression for the Jacobian determinant at the point p of 



shock interaction, (7.7)- (7.8), we get 



det J^(p) = (0oCo - ^0^0) = 0. (7.21) 

Now, the covariant transformation law (2.2), that is, Qap = J^J^g^v, immediately 
implies 



deigaii{p) = {det J^{p)f det g^u{p) = 0. 



This completes the proof of Theorem 7.1 



□ 



We remark that at first there appears to be more than enough freedom to 
choose the free functions $, Q, Z, N of the canonical form to arrange for the 
discontinuous term in the integrability condition to vanish. This together with 



the fact that the derivatives of are uniquely solvable in condition (5.11), lead 



us to believe until the very end that one could construct coordinates in which 
gai3 was C^'^. But at the very last step, taking the limit of the integrability 
constraints to the limit of shock wave interaction p, we find that the crucial 
regularity of $, fl, Z, N, expressing that [^'0/3,7] vanishes at shocks, has the effect 
of freezing out all the freedom in $, Q, Z, N, thereby forcing the determinant of 
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the Jacobian to vanish at p. The answer was not apparent until the very last 
step, and thus we find the result remarkable, and most surprising! 



Shock Wave Interactions are Regularity Singularities in GR; 

THE Full Atlas 



In Theorem 7.1, we proved that C^'^ coordinate transformations of the {t,r)- 
plane, capable of mapping a C°'^ regular SSC metric g^i, over to a metric ga/s in 
C^'^, do not exist. In this section, we extend the result to the full atlas of coor- 
dinate transformations, that is, transformations depending on all four spacetime 
coordinates, thereby proving Theorem 1.1 Recall: 



Theorem 1.1 Suppose thatp is a point of regular shock wave interaction in SSC, 
in the sense of Definition 3.1. for the SSC metric g^^. Then there does not exist 
a C^'^ regular coordinate transformation x°' o [x^)~^, defined on a neighborhood 
M ofp, such that the metric components g^p are functions of the coordinates 
and such that the metric has a nonzero determinant at p. 



The method of proof is as follows, we assume for contradiction the existence 
of coordinates in which the metric is regular, in general the metric does have 



the full ten components and is not of the box diagonal form (2.19). However, in 



spherically symmetric spacetimes one can change to coordinates, such that the 



resulting metric is of box diagonal form (2.19), is regular and is related to 



by a transformation in the (t, r)-plane. This contradicts Theorem 7.1 



For constructing these coordinates we partially follow [25j, (chapter 13). Show- 
ing that this coordinate transformation preserves the metric regularity is our 
own contribution. 



Before we begin with the proof of Theorem 1.1 we need to study the regularity 
of a Killing vector-field of a given metric: 

Lemma 8.1. Let be a Killing vector-field of a metric tensor gij. If in some 
coordinates the metric is in and the Killing vector-field is C^'^ across 
a smooth (timelike) hypersurface S, then, in those coordinates, X^ is also 
regular. 



Proof. Lipschitz continuity across S (c.f. Definition 4.1) implies X., 



sibly discontinuous across S but smooth elsewhere. Denoting with 



across E, continuity of Xj implies (c.f. (4.2)) 



to be pos- 
] the jump 

(8.1) 

for any smooth vector field tangent to S. 

On the other hand, X^ being a Killing vector-field implies that its dual vector, 
Xj, satisfies Killing's equations: 







Xij + Xj^i 



(8.2) 

Taking the jump across S of 
Killing's equations (|8.2[) and using the continuity of the Christoffel symbols F^-, 



where F^- denote the Christoffel symbols of gij. 
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leads to 

[X^A + [X,A = 0. (8.3) 



Define A* = where 5*' denotes Kronecker's symbol. By (8.3), A is 



(point-wise) an antisymmetric matrix. Thus, A is diagonalizeable and has a 



vanishing trace, tr(y4) = 0. Now, by (8.1), A has a three dimensional nullspace 



so that ti{A) = implies all eigenvalues being zero and therefore 

[Xij] = for all z,j e {0,...,3}. (8.4) 



This proves that X* = g^^Xj is in C^, completing the proof of Lemma (8.1 ). □ 



Proof, (of Theorem 1.1) 



Assume for contradiction there exist coordinates , such that the transformed 



metric Qij = ^j^^Q^v is regular and has a non- vanishing determinant at p. In 



general gij is not of the box diagonal form (2.19). We now construct a coordinate 
transformation taking g^j over to a metric of box diagonal form, maintaining the 
metric regularity, in contradiction to Theorem 7.1 



For this, observe that the Killing vectorq^of the SSC metric are smoothp^ and 
are C°'^ regular in any other coordinate system, which can be reached within the 
atlas of C^'^ coordinate transformations. Now, due to the smoothing conditions 



(5.1), a Jacobian lifting the metric regularity to is in C^'^ \ on the shock 
surfaces, while its differentiability elsewhere is negligible as long that it doesn't 
lower the metric regularity. Thus, we assume without loss of generality that 
the Jacobians are smooth off the shock surfaces and consequently the Killing 
vectors are C°'^ across each of the shock surfaces, excluding their intersection. 



Now, since the metric cjij is in C^, Lemma 8.1 immediately shows that those 
Killing vectors are regular on a deleted neighborhood of p, by which we refer 
to a neighborhood of p not containing the intersection of the shock surfaces. 
Given two point-wise linearly independent regular Killing vectors with a 
vanishing Lie-bracket one can choose two vectors, such that those four vectors 
are integrable to coordinate^ where all Killing vectors satisfy (possibly after 
reordering) 

= for j = 0, 1, (8.5) 

while the metric regularity is maintained on a deleted neighborhood. We 
denote these coordinates again with . 

We now construct new coordinates, x^ , such that the resulting metric satisfies 

gal = 0, (8.6) 



"'^'^For ease of notation we subsequently refer to vector fields as vectors. 

"'^^In more detail, there a three smooth generating Killing vectors in SSC and for our purposes 
it suffices to restrict consideration to smooth linear combination of those. 

"'^''A set of vectors is integrable to a coordinate system if and only if their Lie-brackets 
mutually vanish. Now, the condition of vanishing Lie-brackets ties together the regularity of 
the vectors, so that one can always stay within the same class of smoothness. 



38 



M. REINTJES 



for a G {0,1} and / G {2,3}. For ease of notation, we subsequently split up 
indices as a, 6 G {0,1} and k,l,m,n G {2,3}. Now, introduce the coordinates 
implicitly as 

for a G {0, 1}, 



x' 



X 



x\x) for / G {2,3}, 



•7) 



such that X {x) satisfies, (with summation over A; = 2, 3 only). 



dx 



dx"- 



r, 



0, 



with initial data '(xo) for s ome Xq in our coordinate patch. From the 

covariant transformation law, (2.2), we find the metric in the new coordinates, 
gij, satisfying 



9ak 



dx'' ( dx 



dx^ V dx'^ 



-gim + gia 



and using (8.8) we immediately find (8.6) to hold true. Moreover, (8.5) is main- 
tained due to the second equation in (8.7). 



To see that the metric regularity is preserved, write (8.8) in its equivalent 
form 

dx' 



dx° 



where g denotes the 2 — D reciprocal of gik, i.e., g gim = restricting sum- 
mation to I G {2,3}. Since the right hand side in (8.9) is in C^, (2.2) implies 
the metric regularity being maintained. 

We now prove existence of solutions to (8.8), for which it suffices to show that 

(8.10) 



the right hand side of (8.9) allows for commuting derivatives, that is 

{x',x\~x\x)rx\x)) = -^J^ {x\x\~x\x)rx\x)) , 



dx'' 



where we define = g'^'gia- From (8.7), the chain rule and (8.9), we find (8.10) 
to be equivalent to 

KiFl - Kb = F'A - Ka- (8-11) 

We now use Killing's equation to show the validity of (8.11 ). Using (2.8 ), we write 
Killing's equations (8.2) in its equivalent form (summation over a G {0, ...,3}) 

X^'ig^j + X\fg„i + X'^^i,-, = 0. (8.12) 



For indices i = I and j = k, contracting (8.12) with g^^g^^ and using g"^ 



_-in-kmg^^^^ gives 



,l9 



Xk %mn 
9,k 



Xn %lm 
,i9 ■ 



For indices i = k and j = a, contracting (8.12) with g''' gives 

^ ,a — ,kg gma — ^ g gka,m- 



Substituting ( |8.13[ ) into ( |8.14[ ) and recalling = g gia leads to 

,k->- a ~ a,k- 



X'=X^ .Ff - X'T' 



(8.13) 
(8.14) 
(8.15) 
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Taking of (8.15), substituting again (8.14) and using X\^^ — X\^ = we obtain 



+ 



rpm 

b,m a 



Tpk 



a,m b 



b,km a a,km b a,kb 



I ^pl 



b,ka 



a,k b,m 



)(8.16) 



Observe that at each point q and for any antisymmetric 2-form Aik one can find 
a Kilhng vector with X'^(g) = and giiX\{q) = Aik- Therefore, since (8.16) 



holds for all Killing vectors, (8.16) implies (8.11). The existence of a solution 
to (8.9) now follows with standard ODE-methods. (Namely, for fixed I and 
(8.9) is an ODE in x", defining x{x) = x{x"', x^) to be its solution for all x^ with 



the remaining equation in (8.9) (for fixed I). 



h ^ a and for appropriate initial data, (8.10) immediately shows x[x) to satisfy 



In this and the next two paragraphs we show that Qij, in the coordinates x^ 
constructed above, is of box-diagonal form up to a change of (x^, a;^)-coordinates. 
Recall that indices are split as a, 6 G {0, 1} and fc, l,m,n G {2, 3}. Using (8.5) as 



well as (8.6), we find the (a, 6)-components of Killing's equations (8.1) to imply 



X^gab,k 



0. 



Since this holds for all Killing vectors we conclude 

9ab,k = 0, (8.17) 

the (a, 6)-components of the metric are independent of x"^ and x^. 

Denote with cjij the induced metric on one of the spaces of symmetry, which 
we consider fixed and henceforth refer to as the space of symmetry. Then, on 
the spaces of symmetry the metric agrees with the induced metric up to a factor, 
that is, 

gik{x\x\x^,x^) = C{x^,x^)gik{x^,x^), (8.18) 

for some real valued positive function C and k,l G {2, 3}, as we show in the 
remainder of this paragraph. Let q be an arbitrary point lying in our coordinate 
patch. Suppose X'' is a Killing vector satisfying 

X'iq) 



0. 



Now, at q, (8.12) gives 



,i9kj + X^jSki) 



0. 



(8.19) 
(8.20) 



Defining X^^ = {X%^)^^{q) = X\i{q)gkm{q) and 5f = g^^'g^j, is 
equivalent to 



X^M {5\Bf + 5\BT) 



(8.21) 



Thus, since Xmi{q) is antisymmetric due to ( |8.20[ ), the expression in the brackets 
must be symmetric in m and /, that is. 



(8.22) 
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Contracting now with gives 



9 9kj 



Oj 9 9ki, 



which is equivalent to 



9ij — Cgij, 



(8.23) 



for C = \9^^9ki) and since q was arbitrary this holds on the whole coordinate 



patch. Now, C only depends on coordinates x° and x^, since substituting (8.23) 
into (8.12) and using that any arbitrary Killing vector of g is also a Killing 



vector of g, together gives for all 
from which we conclude = for / G {2,3}. This proves (8.18). 



0, 



(8.24) 



Now, (8.6), (8.17) and (8.18) together show that, in the coordinates con- 



structed in the beginning, the metric is given by 

ds"^ = gab{x°, x^) dx"- dx^ + C(x°, x^)gik{x'^, x^) dx^ dx^ , (8.25) 

where summation over a, h runs from to 1, summation over /, k runs from 2 to 
3, C is some function and gij{u) denotes the induced metric on the space of 
symmetry. 

Now, assume without loss of generalitjj^ the space of symmetry has constant 
curvature K = 1, then the induced metric cjijdx^dx^ can be taken over to the 
line element dQ"^ = dd"^ + s\.T?{'d)dLp'^ of the SSC metric we started in by a 
coordinate transformation on the spheres of symmetry alone, (c.f. [25], chapter 



13.2 



Setting t = x^ and r = x^, the resulting metric, 

ds" 



-A{t, r)dP + Bit, r)df'^ + 2E{t, r)dtdf + C{t, r)dQ'^ , (8.26) 

is again in and has the same induced metric on the spheres of symmetry than 
the original SSC metric, g^u- Consequently (8.26) can be taken over to g^i, by 
a transformation in the [t, r)-plane, since by our incoming assumption gij and 
the SSC metric are related by some coordinate transformation. This contradicts 



Theorem 7.1 since the metric (8.26) is in on a deleted neighborhood of p and 
has a non- vanishing determinant at p, which proves of Theorem 1.1 □ 



We close this chapter by stating a straightforward consequence of Theorem 1.1 



asserting that at points of regular shock wave interaction the Einstein equations 
can only hold in the weak sense. 

Corollary 8.2. Suppose that p is a point of regular shock wave interaction in 



SSC, in the sense of Definition 3.1, for the SSC metric g^^. Then there do not 
exist coordinates that can be reached from SSC by a C^'^ coordinate transfor- 
mation, such that the metric is a C^'^ function of the new coordinates, has a 
non-vanishing determinant and solves the Einstein equations point-wise almost 
everywhere in a neighborhood of p. 



"'^'''By assumption of spherical symmetry, the spaces of symmetry have constant curvature, 
so that, without loss of generality, we can choose the space with K = 1. 
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9. The Loss of Locally Inertial Frames 

In this section we prove the non-existence of locally inertial frames around a 
point of regular shock wave interaction. This is in strong contrast to the situation 
of a single shock wave, where such coordinates exist, c.f. [H]. We first clarify 
what we mean by a locally inertial frame. 

Definition 9.1. Let p be a point in a Lorentz manifold and let be a coordinate 
system defined on a neighborhood M of p. We call a locally inertial frame 
around p if the metric Qij in those coordinates satisfies: 

(i) gij{p) = T]ij, where rjij = diag{—l, 1, 1, 1) denotes the Minkowski metric, 

(ii) 9ijAp) =0 for alliJJ G {0,...,3}, 

(Hi) gijki cire in for every compact region contained in M , for all i,j, k,l G 

{o;...,3}. 

We refer to a metric gij, that satisfies {i)-{iii), as a locally Minkowskian {or 
locally flat or locally inertial) metric around p. 

Condition (iii) ensures that, in General Relativity, gravitational effects induce 
second order correction to tensorial equation. That is, the difference of an equa- 
tion in Minkowski space and in a (curved) spacetime is a correction term second 
order in the metric derivatives. (In most literature the metric is assumed to be 
smooth, thus implying condition {iii) of Definition 9.1 ) Now, by Theorem |1.1[ 



there exist distributional and thus unbounded second order derivatives of the 
metric, so that the following Corollary is a straightforward consequence. 

Corollary 9.2. Letp be a point of regular shock wave interaction in SSC in the 
sense of Definition 3.1, then there does not exist a C^'^ coordinate transformation 



such that the resulting metric gij is locally Minkowskian around p. 

Proof. Assume for contradiction that there exist a locally inertial frame around 
the point p of regular shock wave interaction. Then the metric in coordinates 
satisfies gij{p) = rja and thus has a non- vanishing determinant. 



By Theorem 1.1, the metric gij is not in C^. Now, there exist indices i,j,l G 
{0, ..,3} such that gij^i is discontinuous and thus there exist indices i,j,l,k G 
{0, .., 3} such that gij^ik is distributional and therefore not in L°° for any neigh- 
borhood of p. This completes the proof. □ 

Corollary |9 . 2 1 proves that the gravitational metric cannot be locally Minkowskian 
at points of regular shock wave interaction in SSC since unbounded second order 
metric derivatives occur. However, it is presently not clear if there exist coor- 



dinates in which the metric satisfies condition {i)-{ii) of Definition 9.1, or how 
close the metric can be to locally Minkowski. 

10. The Riemann Curvature Tensor is Bounded 

In this section we prove that a regularity singularity is not a naked singularity, 
since the scalar-curvature, the Riemann curvature tensor and all scalar quantities 
incorporating it are bounded. 
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To begin with, recall that the scalar curvature is the trace of the Ricci tensor, 
that is, 

R = R^, 



Taking the trace of the Einstein tensor (2.7) gives 

G% = R% - ^6^,R = -R, 



and taking the trace of the Einstein equations (2.6) we find 



R 



(10.1) 
(10.2) 
(10.3) 



Assuming T^'^ G L°°, (10.3) implies R to be in L°° as well. We conclude that 



the scalar curvature stays bounded at points of regular shock wave interactions 
in SSC. 

However, even though the scalar curvature is bounded, there could be other 
curvature scalars blowing up at points of shock wave interaction, curvature 
scalars of the form 



c = C^-'P-R 



(10.4) 

where R^upa denotes the Riemann curvature tensor and C^^^'^ is some tensor 
with components in L°°, (with respect to some coordinate system). The next 
lemma states that the Riemann curvature tensor and any such scalar is bounded 
for e L°°. 

Lemma 10.1. Assume the energy momentum tensor is hounded almost every- 
where, that is, TP" G Let g^^i, he a regular metric solving the Einstein 
equations (2.6) weakly. Then the Riemann curvature tensor and any curvature 



scalar c, as defined in (10.4), are in L 



Proof. Since CP"P", in (10.4), is in L°°, it suffices to prove that all components 



of the Riemann curvature tensor are in L°°. 

To begin with, we compute the scalar curvature in SSC (using MAPLE): 



R 



Btt ^ AtBt 



+ 



lAr IBr 2 — IB 

+ + (10.5) 



AB ' 2A^B 2A52 rAB rB^ r'^B 

For a Lipschitz continuous metric the first order metric derivatives exist almost 



everywhere, while the second order derivatives in (10.5) have to be taken in a 



weak sense, that is, A„ and B^ are functions in L°° with the defining property 



{A„ - Btt) V 



{Ar^r - Bt^t 



:io.6) 



for all V9 G C^, with denoting the set of smooth functions with compact 
support. By (10.3), we conclude that R G thus (10.5) together with the 



metric being Lipschitz continuous implies 

A„ — Btt 



G 



(10.7) 



where Aj-r and Btt again refer to weak second order derivatives. 

For a metric in C°'^, all components of the Riemann curvature tensor not 
containing second order metric derivatives lie in L°°. In SSC, a straightforward 
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computation shows the only components (up to the symmetries of R^upa) con- 
taining second order metric derivatives to be 

i?oioi = "\ + ' + -^^^^ (10.8) 



where we interpret second order derivatives again in the weak sense. From (10.8) 



and (10.7), we conclude that -Roioi hes in L°°. Finally, since all components of 
R/iupa are in L°° with respect to SSC, the Riemann curvature tensor is in L°° 
in any coordinates that can be reached from SSC within a C^'^ atlas. This 
completes the proof. □ 

By Lemma [10. 1[ even though delta function sources in the second derivatives 
of the gravitational metric tensor exist in all coordinate systems of the C^'^ 
atlas, the Riemann curvature tensor and all curvature scalars remain uniformly 
bounded at regularity singularities due to cancelation. This further clarifies the 
difference between the essential C°'^ singularities in the metric at points of shock 
wave interaction, and the essential singularities at surface layers like the 
"thin shells" introduced in Israel's paper ^U\. Namely, on surface layers, the 
delta function sources in T are the cause of the essential C^'^ singularity in the 
metric g, because second derivatives of g must have distributional sources and 
consequently g cannot be C^'^ in any regular coordinate system. For shock wave 
solutions of G = kT the issue is more delicate because the constraint that G have 
delta function sources is removed and, at first sight, there is no clear obstacle to 
the existence of coordinate systems that smooth the metric to C^'^. In light of 
this, the existence of a non-removable C°'^-singularity at points of shock wave 
interaction is rather subtle and surprising. 

11. Discussion 

In this section we discuss our choice of atlas, the role of spherical symmetry 



for the result of Theorem LI and the possibility of physical effects caused by a 
regularity singularity. 



Theorem Lj^ shows that there is a non-removable lack of regularity within 
the atlas of C^'^ regular coordinate transformations. We believe this proves 
points of regular shock wave interaction to be essential spacetime singularities, 
since we expect that one can not extend the atlas to an even lower regularity and 
achieve a metric-smoothing within the larger atlas. Namely, a Jacobian which 
is Holder but not Lipschitz continuous cannot meet the (necessary) smoothing 



condition (5.4) because it fails to provide enough free parameters [J^^^], said 
differently, a Holder continuous Jacobian does not "mirror" the C^'^ metric reg- 
ularity appropriately. Moreover, we believe that lowering the atlas regularity to 
anything less smooth than would conflict with the weak formulation of the 
Einstein equations, since the metric regularity is too low to shift enough deriva- 
tives onto test functions, due to the quasi-linear structure of the equations. 

Points of regular shock wave interaction are not hidden from observation by an 
event horizon. This opens up the opportunity of direct measurements of effects 
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caused by a regularity singularity. Predicting what such an effect could be is 
an open problem. Supernovae are cosmic events where we believe shock wave 
interactions to take place, namely, the collision of the outer shell of the collapsing 
star with the surface of the (cooled off) core. 

The assumption of spherical symmetry enters our methods in Sections |5] - [7] 
extensively. However, we believe that removing our symmetry assumption does 
not alter our results, since the smoothing condition (5.4) is valid for metrics 
without any symmetries and we expect that (5.4) together with 5.5 are still 
uniquely solvable. But extending our method to the case without any symmetries 
is far from being straightforward. 

Introducing (Navier Stokes type) viscosity terms into the perfect fluid source 
of the Einstein equation causes a smoothing of the discontinuities in the fluid 
variables smearing out the shock profiles. As a result the metric tensor is smooth 
and locally inertial frames exist. However, steep gradients in the fluid sources 
persist and we still expect to recover large second order metric derivatives that 
blow up as the viscosity tends to zero. From this, we conclude that the distribu- 
tional second order metric derivatives could be of physical significance, but the 
loss of locally inertial frames is rather a mathematical curiosity. 



12. Conclusion 

Our results show that points of shock wave interaction give rise to a new kind 
of singularity which is different from the well known singularities of General 
Relativity. The famous examples of singularities are either non-removable singu- 
larities beyond physical spacetime, (for example the center of the Schwarzschild 
and Kerr metrics, and the Big Bang singularity in cosmology where the curva- 
ture cannot be bounded), or else they are removable in the sense that they can 
be transformed to locally inertial points of a regular spacetime under coordi- 
nate transformation, (for example, the apparent singularity at the Schwarzschild 
radius and any apparent singularity at smooth shock surfaces that become regu- 
larized by Israel's Theorem, |10l HH]). In contrast, points of shock wave interac- 
tion are non-removable C"'^ singularities that propagate in physically meaningful 
spacetimes, such that the curvature is uniformly bounded, but the spacetime is 
essentially not locally inertial at the singularity and second order metric deriva- 
tives are distributional. For this reason we call these regularity singularities. 

Since the gravitation metric tensor is not locally inertial at points of shock 
wave interaction, it begs the question as to whether there are general relativistic 
gravitational effects at points of shock wave interaction that cannot be predicted 
from the compressible Euler equations in special relativity alone. Indeed, even if 
there are dissipativity terms, like those of the Navier Stokes Equations]^ which 
regularize the gravitational metric at points of shock wave interaction, our results 
assert that the steep gradients in the derivative of the metric tensor at small 
viscosity cannot be removed uniformly while keeping the metric determinant 

"'^^The issue of how to incorporate a relativistic viscosity that meets the speed of light bound 
is problematic, 
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uniformly bounded away from zero, so one would expect the general relativistic 
effects at points of shock wave interaction to persist. Moreover, points of regular 
shock wave interaction are not hidden from observation by an event horizon, 
which opens up the opportunity of direct measurement of effects that could 
resemble unbounded scalar curvature. We thus wonder whether shock wave 
interactions might provide a physical regime where new general relativistic effects 
might be observed. 



Appendix A. The Integrability Condition 

In this section we review the equivalence of the existence of an integration 
factor and the integrability condition (2.4) on M^. Recall that the Jacobian is 
defined as 



dx^^ 

provided we are given a coordinate transformation from x^ to where the 
indices a = 0, 1 and fi = t,r label a set of four functions. 



Lemma A.l. Let Q be an open set in with coordinates x^ = (t,r). Suppose 
we are given a set of functions J°{x'^) in C'^'^{'D,), away from some curve 
7(t) = {t,x{t)), satisfying det(J°) ^ 0. Then the following is equivalent: 

(i) There exist locally invertible functions x°'{t,r) G C^'^{Q), for a = 0, 1, 
such that g = 

(a) The set of functions J° G C^'^{Q) satisfy the integrability condition 

Proof. The implication from (i) to (ii) is trivial, since (weak) partial derivatives 
commute. We now prove that (ii) implies (i). Without loss of generality, we 
assume that VL is the square region {a,by C M^. For (t,r) G Q, introduce 

a;"(t,r)= / J'^{t,x)dx+ I J^{T,a)dT, (A.3) 

J a J a 

then 

dx^_ 
dr ~ 



follows immediately. Furthermore, using A. 2 we get 

dx°' 



{t,r) = J J^^t{t,x)dx + J^{t,a) 

x{t) pr 

J°,dx+ / Jl,dx + J^{t,a), 

Jx(t) 

where x{t) G (a, r) is the point of discontinuity of J"^. We apply the fundamental 
theorem of calculus to each of the above integrals separately and obtain 

— {t,r) = jnt,r). 
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Moreover, the Inverse Function Theorem yields that the function is bijective 
on some open set, since 




= det ( J^) + 0, 



in Q. □ 
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